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SPECTROSCOPY  AND  STUDIES  OF  ERRORS 

By 
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December,  1976 

Chairman:  James  D.  Winefordner 
Major  Department:  Chemistry 

The  lognormal  distribution  is  a statistics  for  positive 
variables.  It  is  applied  to  the  positive  physical  quantities, 
such  as  transmitted  irradiance,  frequency,  wavelength,  etc. 

It  is  shown  that  the  transmitted  irradiance  to  be  de- 
termined is  not  identical  to  the  quantity  that  is  measured. 

The  way  to  cover  the  error  is  discussed.  The  lognormal 
distribution  explains  the  red-shift,  asymmetry  and  broadening 
of  a spectral  profile.  It  also  explains  why  the  reaction 
rate  constant  is  not  a constant  and  interprets  the  bending 
of  the  analytical  growth  curve. 
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CHAPTER  I 
INTRODUCTION 


lognormal  distribution,  in  its  simplest  form,  may 
be  defined  as  a variable  whose  logarithm  obeys  the  normal  law 
of  probability.1  The  related  literature  is  large,  and  appears 
under  a variety  of  names  - the  Galton-McAlister , Kapteyn  or 
Gibrat  distribution,  the  logarithmico-normal  or  simply  the 
lognormal  distribution.  The  general  properties  of  the  log- 
normal distribution  have  been  well  studied  and  reviewed.1  20 
It  has  found  applications  in  small  particle  statistics,  e.g., 
distribution  of  particle  sizes;  1 50  in  economics,  e.g. , 
distribution  of  prices,  sales  volumes,  incomes,  etc . ; 1 ’ 19 ’ 39-49 
in  biology  and  anthropology,  e.g.,  weights  of  human 
beings;4’6’12’15’ 50-64  in  analysis  of  the  abundances  of 
species,  e.g.,  trace  concentrations  in  rocks;65-^2  in  engi- 
neering, e.g.,  mechanical  or  electrical  lifetimes  of 
18  78—89 

materials;  ’ in  philology,  e.g.,  distributions  of  the 

radical  components  of  Chinese  characters;1,90-94  and  in  physics, 
e.g.,  distributions  of  neutrons  slowed  by  elastic  impacts.95’96 
During  the  past  decade,  the  applications  of  the  lognormal 
distribution  have  included  studies  of  irradiance  distributions 
of  electromagnetic  waves  passing  through  or  scattered  by  in- 
homogeneous media,  and  analysis  of  the  band  shapes  of  molecular 


spectra . 
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The  lognormal  distribution  has  different  forms.1  In  its 
simplest  form,  the  lognormal  is  better  than  the  normal  dis- 
tribution for  positive  quantities.16,18  The  detailed  theory 
of  the  lognormal  distribution  will  not  be  given  here.  However, 
an  example  will  be  given  to  show  why  the  normal  distribution 
is  strictly  applicable  only  if  ...  . negative  measurements 
are  possible."  The  example,  given  here,  is  a study  of  the 
wordlengths  (numbers  of  letters  per  word)  on  page  1 from  a 
book  written  by  H.  A.  Laitinen.130  There  are  304  words  and 
1717  letters  in  that  page.  The  average  wordlength  and  the 
standard  deviation  of  the  wordlengths  are  5.65  and  3.49 
letters  per  word,  respectively.  The  distribution  of  the  word- 
lengths  is  shown  in  Figure  1-1.-  If  it  were  assumed  that  the 
wordlength  followed  the  normal  distribution  with  the  average 
wordlength  being  5.65  and  the  standard  deviation  being  3.49, 
as  shown  in  Figure  1-2,  then  one  would  expect  from  the  normal 
distribution  that  there  should  be  some  words  with  negative 
wordlengths.  However,  negative  wordlengths  are  impossible. 

It  is  clear  in  the  above  example  that  the  normal  distribution 
is  inapplicable  because  negative  measurements  of  wordlengths 
are  impossible. 

In  the  simplest  lognormal  distribution,  according  to  the 
definition,  the  logarithm  of  a positive  variable,  rather  than 
the  variable  itself,  follows  the  normal  distribution.  For 
example,  if  the  pH  values  of  a solution  follow  the  normal  dis- 
tribution, then  the  activities  of  the  hydronium  ions  follow 
the  lognormal  distribution.  One  of  the  fundamental  properties 
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Fig.  1-1  - Actual  Distribution  of  Word  Length  (Mean  = 5.69;  Standard  Deviation  = 3.49) 


aouaaanooo  go  Aouanbaag 


Fig.  1 2 Normal  Distribution  of  Word  Length  (Mean  — 5.69;  Standard  Deviation  = 3.49) 
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of  the  lognormal  distribution  is  that  the  average  of  a variable 
does  not  represent  the  value  of  that  variable  when  there  is 
no  disturbance  of  the  variable.  ’ 3 In  this  case,  for  a 
positive  variable,  the  random  error  will  "not  follow  the  normal 
law  of  error  (the  normal  distribution) , which  is  the  fundamen- 
tal starting  point  for  analysis  of  data.”130  Even  by  taking 
the  average  of  sufficient  observations,  it  is  impossible  to 
overcome  the  scatter  (with  limits)  so  that  the  accuracy  will 
not  be  affected.  To  help  recognize  this  problem  in  the  log- 
normal distribution,  but  not  to  justify  it  theoretically,  an 
example  of  exaggeration  will  be  given  here.  Suppose  there  is 
a solution  with  its  pH  value  being  2.  If  the  measured  pH 
values,  due  to  some  unknown  factors,  are  1,  2,  and  3,  the 
average  pH  value  can  be  reported  as  2,  and  the  average  activity 
of  the  hydronium  ions  may  be  reported  as  0.01.  However,  from 
the  same  pH  values,  1,  2,  and  3,  the  activities  of  the  hydro- 
nium ions  can  be  expressed  as  0.1,  0.01,  and  0.001,  respectively. 
Therefore,  the  average  activity  of  the  hydronium  ions  may  also 
be  reported  as  (0.1  + 0.01  + 0.001)/3  = 0.037,  which  is  quite 
different  from  0.01  (pH  = 2).  The  same  problem  might  happen 
in  absorption  spectroscopy,  due  to  the  linear  relationship 
between  the  concentration  of  the  absorbing  sample  and  the 
logarithm  of  the  transmitted  intensity  (or  photon  flux).  If 
the  sample  concentration  in  the  optical  region  of  the  sample 
cell  approximately  follows  the  normal  distribution,  then  the 
transmitted  intensity  would  follow  the  lognormal 

q 7 _ i 99 

distribution . 
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The  objective  of  Part  I of  this  thesis  is  to  apply  the 
lognormal  distribution  to  analytical  optical  spectroscopy  and 
to  the  study  of  indeterminate  errors  in  analytical  optical 
spectroscopy.  Part  II  covers  the  determinate  errors  caused 
by  the  multiple  reflections  within  the  cell. 


PART  I 


APPLICATIONS  OF  THE  LOGNORMAL  DISTRIBUTION 
TO  OPTICAL  SPECTROSCOPY 
AND  STUDIES  OF  INDETERMINATE  ERRORS 


CHAPTER  II 


LOGNORMAL  DISTRIBUTION  OF  THE  IRRADIANCES  OF  AN 
OPTICAL  WAVE  AFTER  PASSING  THROUGH  A RANDOM  MEDIUM 


The  problem  will  be  first  treated  with  a binomial  model 
and  then  with  a random  walk  model.  Eventually,  the  lognormal 
distribution  of  the  transmitted  irradiances  (intensities  or 
photon  fluxes)  of  an  optical  wave  passing  through  a turbulent 
medium  is  derived  without  the  use  of  the  complicated  theory 
of  electromagnetic  waves,  assuming  the  numbers  of  the  filters 
(absorbers  or  optical  attenuators)  in  the  optical  region 
follow  the  normal  distribution.97'122  Expressions  for  in- 
accuracies, signals  and  noises  are  derived. 


Binomial  Model 

In  order  to  understand  the  basic  properties  of  the  log- 
normal distribution,  it  is  helpful  to  picture  an  optical  medium 
as  composed  of  four  adjoining  cells  for  the  illustration. 

Under  the  condition  of  no  turbulence,  each  of  the  four  cells 
is  assumed  to  contain  a filter  (or  absorber) , with  the  same 
transmittance,  t.  The  light  source,  S,  assumed  to  be  mono- 
chromatic, gives  off  a parallel  beam  of  radiation  with  an 
incident  irradiance  IQo  (s  1) . After  passing  through  the 
homogeneous  medium  with  one  filter  in  each  of  the  four  cells 
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Hypothetical  Homogeneous  Medium 
S:  Light  Source;  D:  Detector; 
l-oo:  Incident  Irradiance;  IQ  Transmitted 
Irradiane;  t:  Transmittance  of  the  Filter 


9 


as 


shown  in  Figure  2-1,  the  irradiance  of  the  transmitted 

_ T , -lx  . 


light,  IQ  (s  ),  is  given  by 
I 


■o  - t4 


(2-1) 


The  absorbance  of  the  homogeneous  medium,  A , is  given  by 

I 

OOn  , , 

(2-2) 


Ao  “ log  = -4  log  t 

o 


However,  under  the  condition  of  turbulence,  each  cell  is 
assumed  to  be  statistically  identical  and  independent,  and  to 
have  equal  probability  of  either  gaining  an  extra  filter  or 
losing  one.  The  combinations  of  the  possible  occurrence  of 
gaining  or  losing  a filter  in  each  of  the  four  cells  are  shown 
in  Table  2-1,  where"+"denotes  that  the  cell  gains  one  more 
filter  or  contains  two  filters,  and"-"denotes  that  the  cell 
loses  one  filter  or  contains  no  filter.  In  Table  2-2,  the 
relative  frequencies  of  occurrences  and  the  expressions  for 


the  transmitted  irradiances  in  the  turbulent  medium  are  given. 
The  relative  frequency  versus  the  overall  transmittance  of 
the  turbulent  medium  in  Table  2-2  is  plotted  in  Figure  2-2 
f°r  t - 0.7  and  IQO  = 1 , as  an  example.  The  average  trans- 
mitted irradiance,  f ave . of  the  turbulent  medium,  is  given  by 


ave 


(t8  + 4 t6  + 6 t4  + 4 t2  + 1) 


= (1  + t2)4 

TB  ^oo 


(2-3) 


So  the  absorbance  of  the  inhomogeneous  medium,  A.  , is  given  by 
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Table  2-1 

Possible  Combinations  and  Number  of 
Filters  in  the  Turbulent  Medium 


Possible  Combination 


Number  of  Filters 
in  the  Medium 


+ 

+ 

+ 

+ 

+ 

+ 

+ 
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+ 

+ 

+ 
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+ 
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+ 

+ 

+ 

+ 


+ 

+ 
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+ 

+ 


+ 
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+ 


+ 


+ 


+ 


+ 


8 

6 

6 

4 

6 

4 

4 

2 

6 

4 

4 

2 

4 

2 
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Table  2-2 

The  Relative  Frequency  of 
Occurrence  and  the  Transmittance 


Relative  Frequency  Transmittance  Absorbance 


WW 

l0s  dec/1 

1 

t8 

-8  log  t 

4 

t6 

-6  log  t 

6 

t4 

-4  log  t 

4 

t2 

-2  log  t 

1 

1 

0 
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Fig.  2-2  - Relative  Frequency  and  Transmittance  of  a Turbulent  Medium. (t 
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in 


- log  (jS  . -4  log  (i+-d) 

ave 


For  t - 0.7,  we  have,  from  Equations  (2-1),  (2-2), 


and  (2-4), 
I 

ave 


A. 


in 


I 

o 


A 


o 


0.31 


I 

oo 


0.51 


0.24  I 

oo 

0.62 


(2-4) 


(2-3), 


It  is  clear  that  the  average  transmitted  irradiance  I 

’ ave 

(s  ) , of  an  inhomogeneous  medium,  is  greater  than  the  trans- 
mitted irradiance  of  a homogeneous  medium,  I and  the  ab- 
sorbance of  an  inhomogeneous  medium,  A^,  is  smaller  than  the 
absorbance  of  a homogeneous  medium,  A . Therefore,  in  the 
spectroscopic  determination  of  the  concentration  such  as  in 
a turbulent  flame  or  in  a stop-flow  technique,  the  turbulence 
will  result  in  a lower  value.  The  relative  error  of  the  ab- 
sorbance, E^,  is 


log  (1  +t2)  - log  (2t) 
log  t 


(2-5) 


The  geometrical  mean  of  the  transmitted  irradiance,  I is 

g 


I = [t8  x (t6)4  x (t4)6  x (t2)4  x 1]1/16  I 

6 OO 

= t4  I = I 
OO  o 


(2-6) 


Thus  the  geometrical  mean  of  the  transmitted  irradiances  in 
turbulent  medium  is  equal  to  the  transmitted  irradiance  in  the 
homogeneous  medium.  The  relative  frequency  versus  log  (I/I  ) 
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is  plotted  in  Figure  2-3.  It  can  be  seen  that  in  Figure  2-3 
the  distribution  is  symmetric  about  the  mean,  while  the  dis- 
tribution in  Figure  2-2  is  not  symmetric.  From  the  mean  in 
Figure  2-3,  the  absorbance,  ^ given  by 

Aave  = "TF[log  (t8)  + 4 loS  (t6)  + 6 log  (t4)  + 4 log  (t2)] 
= -4  log  t = Ao  (2-7) 

Thus  the  average  absorbance  is  equal  to  the  absorbance  of  the 
homogeneous  medium.  Equation  (2-6)  or  (2-7)  should  be  used 
to  avoid  errors  due  to  the  turbulence  of  the  medium.  In  prac- 
tice, Equation  (2-6)  is  inconvenient.  It  is  advisable,  from 
Equation  (2-7),  to  take  the  average  after  the  signals  are  con- 
verted to  absorbances.  The  standard  deviation  of  the  trans- 
mitted irradiances,  a^,  is 


a- 


(1  + t4)4 


(1  + t2)8  % 

256  oo 


(2-8) 


'I  " L T6 

Similarly,  if  the  turbulent  medium  consists  of  n cells, 
the  average  transmitted  irradiance,  I . and  the  standard 

cl  v 6 

deviation  of  the  transmitted  irradiances,  Oj. , are  given  by, 
respectively 


'ave 


(1  + t2)n 
0n 


'oo 


(2-9) 


and 


°i  = 1 


(1  + c4)n 


-.n 


22n  oo 


The  absorbance  of  the  homogeneous  medium,  AQ , is 


(2-10) 


tO 
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00 
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Fig.  2-3  - Distribution  of  the  Transmitted  Irradiances  on  the  Logarithmic  Scale 
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Ao  * -n  loS  t (2-11) 

The  absorbance  of  the  inhomogeneous  medium,  A.  is 

in  ’ 

Ain  = loS  = n lo§  (— 2 ?)  (2-12) 

ave  1 + 

Ao  and  Ain  versus  n for  t = 0.7  are  plotted  in  Figure  2-4. 

The  relative  error  of  the  absorbance,  E^,  is  given  by 

E.  = -ill- jo  = log  (1  + t2)  - log  (2t)  (2-13) 

A A0  rsg-t  — 

It  can  be  seen  that  the  relative  error  of  the  absorbance  is 
independent  of  the  numbers  of  the  cells  in  the  medium,  n,  and 
that  it  decreases  as  the  transmittance  of  each  filter,  t, 
becomes  larger  as  shown  in  Figure  2-5.  If  the  optical  activity 
coefficient  of  the  filter  in  the  turbulent  medium,  a,  can  be 
defined  as 


a 


m 


o 


then 


a = l°g  (1  + t ) - log  2 
log  t 


(2-14) 


It  is  clear  that  the  activity  coefficient  is  dependent  on  the 
transmittance  of  each  filter.  The  ratio  of  the  average  trans- 
mitted irradiance  of  the  inhomogeneous  medium,  I , to  that 


of  the  homogeneous  medium,  Iq , is  given  by 
Iave  _ (1  + t2)n 


o 


2n  tn 


(2-15) 


^ave^o  versus  the  transmittance  of  each  filter,  t for  various 


Number  of  Cells 

Absorbances  of  a Homogeneous  and  an  Inhomogeneous  Medium  (t 


Relative  Error  of  Absorbance 


Fig.  2-5  - Relative  Error  of  Absorbance 
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media  is  plotted  in  Figure  2-6.  It  can  be  seen  that  the 
irradiance  ratio  increases  as  the  numbers  of  cells  become 
larger,  or  as  the  transmittance  of  each  cell  becomes  smaller. 
The  signal- to -noise  ratio  of  irradiance,  S/N,  is  given  by 


S 

N 


I 


ave 


(1  + t2)n 

t(l  + t4)n2n  - (1  + t2)2n]^ 


(2-16) 


S/N  versus  n,  for  various  t are  plotted  in  Figure  2-7,  from 

which  we  can  see  that  the  signal-to-noise  ratio  decreases 

as  n increases  or  t becomes  smaller.  For  n = 10,  t = 0.9, 

the  relative  frequencies  of  occurrence  versus  I/I  and 

oo 

l°g  (I/I00)  are  plotted  in  Figures  2-8  and  2-9,  respectively. 
The  distribution  of  I,  in  Figure  2-8.  is  skewed  to  the  right, 
in  other  words,  positively  skewed,  or  the  dominant  tail  of 
the  distribution  is  on  the  right  while  the  distribution  of 
log  (l/Io0)  as  shown  in  Figure  2-9,  is  symmetric  about  the 


mean. 


Substituting  t = 0.7  into  Equation  (2-13)  gives 

EA  = —a  A°  = -0-18 

o 

The  above  equation  indicates  that  in  absorption  spectroscopy, 
the  absorbance  would  have  a lower  value  due  to  the  turbulence 
of  the  optical  medium.  Substituting  n = 4 and  t = 0.7  into 
Equation  (2-15)  gives 

Iave  - OJ  I„ 


The  above  equation  indicates  that  in  fluorescence  spectroscopy, 
the  fluorescence  signal  would  have  a lower  value  because  more 
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Fig.  2-7  - Dependence  of  the  Signal-To-Noise  Ratio  on  the 
Number  of  Cells  and  the  Transmittance  of  the 
Filter 
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Fig.  2-9  - Distribution  of  the  Logarithm  of  the  Transmitted  Irradiance  in  the 
Turbulent  Medium  (n  = 10;  t = 0.9;  Dotted  Curve  is  the  Envelope  of 
the  Distribution) 
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light  would  pass  through  a turbulent  medium  than  a homogeneous 
medium. 

Now,  let  us  consider  the  inhomogeneous  medium  as  con- 
sisting of  (n  + h)  cells  of  which  n cells  suffer  turbulence, 
while  h cells  suffer  no  turbulence,  then 


I 

ave 

_ (1  + t2)n  rh  T 

2n  oo 

(2-17) 

aI  = 

rd  + tV  (1  + t2)2n  .%  rh  T 

2n  22n  oo 

(2-18) 

I 

o 

tn  + h i 

OO 

(2-19) 

'''ave 

Xo 

_ (1  + t2)n 

2ntn 

(2-20) 

I 

g 

Xo 

(2-21) 

S 

(1  + t2)n 

(2-22) 

N 

[(1  + t4)n2n  - (1  + t2)2n]% 

> 

o 

ii 

- (n  + h)  log  t 

(2-23) 

Ain  = 

-h  log  t - n log  (X  2 ) 

(2-24) 

and 


" Ao  = log  (1  + t2)  - log  (2t) 
Ao  (n  + h)  log  t 


(2-25) 


It  is  interesting  to  note  that  Equations  (2-15)  and  (2-20)  are 
identical,  and  that  Equations  (2-16)  and  (2-22)  are  identical, 
even  though  they  apply  to  different  conditions.  We  can  see 
that  Iave/Io  and  S/N  are  independent  of  the  numbers  of  cells, 
h,  which  suffer  no  turbulence.  The  geometrical  mean  of  the 
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transmitted  irradiances,  I is  equal  to  I . The  relative 

o O 

error  of  absorbance,  E^,  decreases  as  t increases. 

However,  the  binomial  model  has  the  mathematical  in- 
convenience in  computations.  Because  the  experimental 
evidences  indicate  a lognormal  distribution  of  irradiances, 
an  approximate  approach  to  the  lognormal  distribution  will 
be  given  in  the  following  section  based  on  the  random  walk 
model . 


Random  Walk  Model 


We  will  assume  the  medium  as  composed  of  k cells,  of 
which  n cells  suffer  turbulence . and  the  rest  (k  - n)  of  the 
cells  suffer  no  turbulence.  To  start  the  discussion,  con- 
sider that  no  cells  suffer  turbulence,  so  the  transmitted 
irradiance,  I is  given  by 

lo  * tk  I00  (2-26) 

If  only  one  cell  suffers  turbulence,  the  transmitted  irradiance, 
I-p  is  given  by 

k ± 1 „ 


- c' 


oo 


= t±al  I 


where  a-^  = ±1. 


If  j cells  of  the  total  k cells  suffer  turbulence,  then 

x _ = tk±l±l....±l±l 

j ioo 
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j 

-£  a. 

= t 1 J I 


where  aj  = ±1. 

In  our  model,  n cells  of  the  total  k cells  are  assumed 
to  suffer  turbulence,  so  the  irradiance,  I,  reaching  the  sur- 
face of  the  detector  is  given  by 


n 

-T.  a. 

I = t 1 J I 


or 


n 

E a. 
i J 


o 


1 o 

ThTT  ln  I 


(2-27) 


(2-28) 


For  simplicity,  define 


n 

N = Z a. 
i J 


(2-29) 


where  N is  the  summation  of  n random  numbers  of  ±1,  then 

N = nrr  ln  -r  (2-30) 

The  combinations  of  n random  numbers  of  ±1  will  be  treated 
with  the  random  walk  model  as  shown  in  Figure  2-10,  using  the 

1 Ol 

Galton  board  for  illustration.  In  this  apparatus,  small 

balls  of  equal  diameter  run  down  a slightly  inclined  plane, 
passing  through  n horizontal  rows  of  pins  arranged  in  such  a 
way  that  a network  of  equilateral  triangles  is  formed.  All 
balls  are  released  at  the  same  spot,  x = 0,  above  the  central 
pin  of  the  first  row.  The  diameter  of  the  balls  is  slightly 
smaller  than  the  distance  between  two  neighboring  pins,  the 
latter  distance  may  be  taken  as  two  units.  It  is  assumed  that 
each  ball  when  hitting  a nail  (of  any  row)  will  deviate  by 
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Fig.  2-10  - The  Galton  Board 
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one  unit  to  the  right  or  to  the  left  with  equal  probabilities, 
1/2.  Each  ball,  after  passing  through  the  n rows  of  pins, 
will  come  to  rest  in  one  of  the  compartments  (the  width  of 
each  being  equal  to  two  units)  at  the  bottom  of  the  board. 

The  probability,  P2(N),  that  a ball  arrives  in  the  compartment 
the  acscissa  of  whose  midpoint,  N,  which  is  the  summation 
of  n random  numbers  of  ±1,  is  given  by  131,132 


P2(N)  = l%(n  + N)]?"[%(n  - N)  ] ! ^ (2-31) 

where  N,  the  summation  of  n random  numbers  of  ±1,  is  a 

positive  or  negative  number  equal  to  the  following131,132 

-n,  -n  + 2,  -n  + 4,  n - 4,  n - 2,  n 

When  N is  sufficiently  large,  Equation  (2-31)  becomes  ap- 
131  112 

proximately  ’ 

2 

p2(N)  = JL-  exp  (-^)  (2-32) 

A numerical  comparison  of  the  two  Equations  (2-31)  and 
(2-32)  is  duplicated  in  Table  2-3  for  n = 10. 131  The  prob- 
abilities , calculated  from  Equation  (2-31),  are  shown  by  the 
histogram  in  Figure  2-11.  When  the  width  of  each  compartment 
is  unity,  we  have  the  probability132 


P(N)  = 


p2(n> 

C o p2(N)  dN 


/2frn 


( N N 

exp  (-2^) 


From  Equations (2-30)  and  (2-31),  we  have  the  probability 


(2-33) 


Table  2-3 


The  Problem  of  Random  Walk  for 


Position 

N 

0 

2 

4 

6 

8 


Probability 
from  Eg.  (2-31) 
0.24609 
0.20508 
0.11715 
0.04374 
0.00977 


n = 10 

Probability 
from  Eg.  (2-32) 
0.252 
0.207 
0.113 
0.042 
0.010 


10 


0.00098 


0.002 
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Fig.  2-11  - The  Probability  in  Random  Walk 
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density  function  of  the  transmitted  irradiances133,13^ 


P(I)  = 


Empirically, 


P(N)  |g| 

1 

/2-jrn  In  t I 

the  lognormal 


exp  [ ^ — In2  (I  (2-34) 

2n  lnz  t o 

density  function  takes  the  form 


P(I)  = — exp  — K In2  (f-)  (2-35) 

/2tt  0 I Io 

where  0 is  the  characteristic  parameter  of  the  inhomogeneous 
123  — 

medium.  Let  0 = -/n  In  t,  and  Equation  (2-35)  can  be 

derived  from  Equation  (2-34) . 

Equation  (2-35)  is  plotted  in  Figure  2-12  for  various 
values  of  3,  and  in  Figure  2-13  for  various  values  of  I 
The  distributions  of  irradiances  in  both  figures  are  skewed 
to  the  right.  However,  the  probability  density  function  of 
the  distributions  of  the  logarithms  of  irradiances,  given  by 

p(m  i)  - p(i)  ijn^ri 


1 2 T 

= — = — exp  [ j In  (y-)] 

/2ir  3 23  o 

is  symmetric  about  the  mean  and  equal  to  In  I 

1 R 

(2-35) , the  rth  moment  of  I is  given  by 


(2-36) 
From  Equation 


Ufl)  = r Ir  P(I)  dl 


2 2 

= exp  (^-)  (2-37) 

Hence,  the  arithmetic  mean  or  the  expected  value  of  I,  I 

r ’ ave ’ 


P(I) 
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Fig.  2-12  - Probability  Density  Function  for  Various  Values 
of  I (Io  = 1) 


P(I) 
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I 


Fig.  2-13  - Probability  Density  Function  for  Various  Values 
of  IQ  (6  = 0.7) 
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is 

62 

^ave  - Xo  exp  (2-38) 

The  above  equation  indicates  that  the  average  (measured) 

transmitted  irradiance  of  an  inhomogeneous  medium  I is 

’ ave  ’ 

greater  than  the  transmitted  irradiance  of  a homogeneous 
medium,  Iq . To  overcome  the  error,  due  to  the  inhomogeneity 
of  the  medium,  g should  be  measured  and  this  requires  addi- 
tional work.  From  Equation  (2-36),  the  average  value  of 
the  logarithms  of  the  transmitted  irradiances  of  an  inhomo- 
geneous medium,  (In  i)ave>  is  given  by 

I = oo 

(ln  I)ave  = J ln  1 P(ln  I)  d(ln  I)  = In  I 
1=0  ° 

Rearranging  the  equation  above  gives 

1^  = exp  [(In  I)  ] 

o r l , aveJ 

The  above  equation  shows  the  error  can  be  avoided  if  one 
takes  the  average  value  of  the  logarithms  of  the  transmitted 
irradiances  of  an  inhomogeneous  medium. 

The  standard  deviation  of  intensities,  a^,  is  given  by18 

°i  = - W2  pw  dl 

2 

= 1Q  exP  (f~)  [exp  (g2)  - 1]^ 

= Iave  Cexp  “ U ^ (2-39) 

It  is  interesting  to  note  that  the  standard  deviation 
of  the  intensity  is  proportional  to  the  average  intensity. 


(2-40) 


The  first  shape  factor,  y^,  or  skewness,  given  by 


/„  (I  - Xave>  fd)  dI 


+ 


> 0 


2 > 

where  y = [exp  (8  ) - 1]  2 is  for  the  simplicity  of  the 
expression,  indicates  that  all  the  lognormal  distribution  is 
skewed  to  the  right  as  shown  in  Figures  2-12  and  2-13. 18  The 
second  shape  factor,  or  kurtosis  (peakedness),  given  by 


y2 


_ oo 

\ 4 (I  - Lve)4  pw dI 

°I 


- 6y6  + 15y^  + 16y2  + 3 > 3 (2-41) 

shows  the  lognormal  density  function  is  more  peaked  than  the 

S 18 

normal  density  function.  ’ From  Equations (2-40)  and  (2-41), 
it  is  seen  that  for  small  values  of  8,  the  shape  factors  of 
the  lognormal  function  are  close  to  their  normal  values,  0 
and  3,  respectively.  The  approximate  normal  distribution, 
P^(I).  obtained  from  the  lognormal  distribution,  when  8 is 
small,  is  given  by 

2 

pNU)  = — -1 exp  [-L1  ^n)  ] (2-42) 

/2tt  8 I 28Z  I2 

o o 

It  will  be  shown  that  Equation  (2-42)  is  an  approximation 
of  Equation  (2-35),  when  8 is  small.  Let  the  difference 
between  Equations  (2-35)  and  (2-42)  be 
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D (I)  = 


~~ exp  [ j-  In2  (f-.) 

/2tt  g I 2gz  io 


1 r 1 
— = exp  [ ~ y (I 

/2ir  3 I 2gz  Iz 

o o 


I )2] 

o'  J 


(i)  When  1=1, 


D (I)  = 0 


(ii)  When  I is  in  the  vicinity  of  I , 

J o 


I = I 


'o 


(2-43) 


(2-44) 


so 


2,1 


I - I 


ln^  (y— ) = ln^  (1  + — = — -)  = 


(i  - I0) 


o 


(2-45) 


Combining  Equations  (2-43)  , (2-44)  and  (2-45)  , gives 
D (I)  = 0 

(iii)  When  I is  far  from  I , both  the  exponential  terms  in 
Equation  (2-43)  fall  to  zero,  and  so 
D(I)  = 0 

From  Equation  (2-35),  the  most  probable  irradiance, 

I „ . obtained  from 
max 


dP(I)  = n 

dl  U 


(2-46) 


is  given  by 


Wax  ' W exP  <-e2>  (2-«) 

From  Equations  (2-38)  and  (2-47) , it  can  be  seen  that  the 
average  irradiance  shifts  towards  right  and  the  most  probable 
irradiance  shifts  towards  left  as  g increases.  The  maximum 
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of  the  lognormal  density  function,  P , defined  as 

max 

pm„v  = p(i  > 

max  v max' 

is  given  by 

Pmax  " 7F"1  T exP  O2)  (2-48) 

v Zu  3 IQ 

The  above  equation  is  plotted  in  Figure  2-14.  It  can  be 
seen  that  Pmax  decreases  and  then  increases  as  3 increases. 

The  signal-to-noise  ratio,  S/N,  is  given  by 

I = = [exP  (2-49) 

Dependence  of  Signal-to-Noise  Ratio 
on  the  Observation  Time 


We  will  take  a look  at  a simple  example,  before  going 
into  a general  discussion.  Let  us  consider  the  throwing  of 
a perfect  hexahedral  die.  Each  throwing  is  independent,  so 
the  probabilities  of  the  numbers  of  pips  are  equal  to  1/6 
as  shown  in  Figure  2-15  (a).  If  two  identical  dice  are 
thrown  at  a time,  then  the  probabilities  of  the  average 
number  of  pips  are  shown  in  Figure  2-15  (b).  It  can  be  seen 
that,  in  the  case  of  throwing  two  dice,  the  probabilities  of 
the  average  number  of  pips  are  more  centrally  distributed 
about  the  mean  than  those  when  one  die  is  thrown  at  a time. 

In  the  analogy  to  the  throwing  of  two  dice,  let  us 
consider  the  average  of  the  positions  of  T balls,  when  T 
balls  are  released  in  the  Galton  board.  The  variance  of  the 
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Fig.  2-14  - The  Maximum  Probability  Density 


Probability  Probability 
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2 3 

Number  of  Pips 

(a) 


4 


Fig.  2-15  - The  Probability  of  (a)  Throwing  one  Die 

(b)  Throwing  two  Dice 
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average  of  the  positions  of  T balls,  n^,  is 


nT  = ^ (2-50) 

Therefore,  the  probability,  PT(N) , is  approximately  given  by 

2 

Pt(N)  = exp  (-^-) 


/2 


Trnn 


/2 


/T  . TNZ. 

exP 


frn 


(2-51) 


If  the  observation  time  of  the  detector  is  increased  T time, 
the  probability  density  function,  PT(I) , is  approximately 
given  by 


PT(I)  = 


pt(n>  igfi 


/T  , T , 2 , I * .. 

— = exp  [ 7 In  (= — )] 

/2tt  $ i 2g  oT 


(2-52) 


where  IQrp  is  the  geometric  mean  of  the  transmitted  irradiances 
The  signal- to-noise  ratio,  S/N,  when  the  observation  time  is 
increased  T times,  is 


S 

N 


= [exp  (^p)  - 1] 


% 


(2-53) 


Equation  (2-53)  is  plotted  in  Figure  2-16.  When  T is 
sufficiently  large,  then 


S 

N 


/T 

B 


(2-54) 


From  Equations  (2-40) , (2-41)  and  (2-52)  the  probability 
density  function,  P^(I) , approaches  a normal  function  as  T 
increases,  because  g//T  becomes  smaller. 


Fig.  2-16  - Dependence  of  Signal-To-Noise  Ratio  on  the 
Observation  Time  (3=1) 
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As  the  observation  time  increases  T times,  the  average 
irradiance  should  remain  the  same,  i.e., 


or 


then 


J I P (I)  dl  = / I P(I)  dl 


IoT  exp  (7t)  = Tn  exP  (— ) 


■oT 


-2T 


o 


(2-55) 


(2-56) 


Io  exp  (2 It) 


The  most  probable  irradiance,  I 
time  increases,  derived  from 


T , max 


(2-57) 

, when  the  observation 


dPT(I) 

dl 


is  given  by 


= 0 


(2-58) 


"T,  max 


= IqT  exp  (-^-) 

2 2 

= I exp  C^-  - ^-) 
o P 4 2T  ' 


(2-59) 


Dependence  of  the  Signal-to-Noise  Ratio 
on  the  Aperture 


Similarly,  when  the  observation  time  is  increased  T 
times  and  the  aperture  is  increased  W times,  and  if  the  source 
is  homogeneous  and  bigger  than  the  aperture,  then  the  proba- 
bility density  function  of  the  transmitted  irradiances 

/TW  r TW  i 2,1  N , 

- exp  [ 7 In  (= )] 

23  oTW 


PTW<^I^) 


/2tt  3 I 


(2-60) 


43 


where  IqTW  is  the  geometrical  mean  of  the  transmitted  irra- 
diances  given  by 

2 2 

IoTW  Io  ~ TTW^  (2-61) 

The  signal-to-noise  ratio  is 


S 

N 


[exp  (^)  - 1] 


-% 


When  T or  W is  sufficiently  large,  we  have 


(2-62) 


S 

N 


/TW 

3 


(2-63) 


From  the  above  equation,  we  see  that  the  signal-to-noise 
ratio  is  proportional  to  the  square  root  of  the  observation 
time  or  the  aperture. 


CHAPTER  III 

THE  LOGNORMAL  APPROACH  TO  SHOT  NOISE 

In  the  previous  chapter,  only  the  noise  caused  by  the 
randomness  of  the  medium  was  considered.  In  this  chapter, 
the  noise  due  to  the  nature  of  photons  caused  by  the  light 
passing  through  a homogeneous  medium  is  considered. 

Bernoulli's  Distribution133,13^ 


The  consideration  here  is  that  the  medium  is  homogeneous 
with  a transmittance,  t^,  and  that  a photon  counter  is  employed 
with  a quantum  efficiency,  n • Therefore,  the  probability  of 
a single  photon,  p,  passing  through  the  medium  and  being 
detected  is  given  by 

p = V1  (3-D 


The  probability  of  x photons  being  detected  by  the  photon 
counter,  Pm(x) , when  m photons  are  emitted  from  a stable  light 
source  and  are  travelling  toward  the  detector  through  the 
homogeneous  medium  is  given  by 


- rum  - x),  Px  (1  - p)"'x 

and  the  expected  (or  average)  number  of  photons  being 
by  the  detector,  E(x),  is133 


(3-2) 

detected 
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E(x)  = m p (3-3) 

Also,  the  standard  deviation  of  the  numbers  of  the  detected 
photons  is 

°x  = /mp(l  - p)  (3_4) 

When  m is  sufficiently  large,  it  can  be  shown  that  Equation 

I 04 

(3-2)  becomes  asymptotically 


Pm(x)  = 
m ' 


1 


r (x  - mp;  , 

exP  [-2i-pTi  - p)  1 


(3-5) 


/2?  /mp  (1  - p ) 

If  x is  treated  as  a continuous  variable,  we  have  approximately 
E (x)  = mp  (3-6) 

and 


a. 


x 


= /mp  (1  - p) 


(3-7) 


The  approximation  is  pretty  good  for  small  m (_>  5)  , as  long 
as  p is  not  too  far  from  1/2.  When  m is  very  large  and  p is 
very  small,  Equation  (3-2)  approaches  the  Poisson  distribution 


given  by 


137-139 


Pm(x)  = "j/? ) exP  (~mP) 

with  mean 


(3-8) 


E (x)  = mp  (3-9) 

and  standard  deviation 

a = /mp  (3-10) 

When  n is  extremely  large,  Equation  (3-8)  approaches  the 
normal  distribution 


pm(x>  = 


m 


/Zu/mp 


exp  [ 


(x  - mp) 
2mp 


(3-11) 


137  138 

and  the  mean  and  the  standard  deviation  remain  the  same.  ’ 
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Bernoulli's  distribution,  given  by  Equation  (3-2),  is 
very  tedious  for  calculation  and  manipulation.  The  normal 
distribution,  given  by  Equation  (3-5),  is  very  convenient 
for  calculation,  but  for  small  p it  does  not  have  the  asymmetry 
of  distribution.  Poissonian  distribution,  given  by  Equation 
(3-8) , has  the  asymmetry  of  distribution  and  is  less  tedious 
than  Bernoulli  s distribution,  but  it  is  still  inconvenient 
for  calculation  and  manipulation,  because  of  the  factorial  term 
in  it.  An  approximate  expression  for  Bernoulli's  distribution 
will  be  given  for  having  both  the  asymmetry  and  the  convenience 
of  calculation. 

Lognormal  Approximation  to  Bernoulli's  Distribution 

When  m is  very  large,  Stirling's  formula  is  given  by139 
In  m!  = m In  m - m + In  /2frm  (3-12) 

Equation  (3-2)  and  Stirling's  formula,  when  both  m and  x are 
very  large,  become 


ln  p U) 

m 


(3-13) 


Let  us  define  a parameter,  u given  by 


x 

mp 


(3-14) 


where  q = 1 - p . 
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When  m is  very  large,  Equation  (3-14)  becomes 


5-  = 1 + 

mp  ■ mp 


and  we  have 


In  — = In  (1  + u^) 

mp  • mp 


= u 


mp 


2 

u_jq 

2mp 


From  Equations  (3-15)  and  (3-16), 


X ln  m^  = u,/^Pq  - \ u2q 

From  Equation  (3-15) , 

m - x = mq  - u/mpq 
so , 

In  (1  - 

- UvII  - UJL 

mq  2mq 


In 


m - x 
mq 


then, 


(3-15) 


(3-16) 


(3-17) 


(3-18) 


(3-19) 


(m  - x)  In  ^ x = - Uv4npq  - ^ u2p  (3-20) 

Substituting  Equations  (3-17)  and  (3-20)  into  Equation 
(3-13)  yields 


In  Pm(x) 
Therefore , 


2 

u 


2 


+ 


m 

2ttx  (m 


30 


(3-21) 


Pm(x) 

m 


r m 
2iTx(m  - x) 


2 

exp  ( - ^ — ) 


1 / mx  1 
x’zit  (m  - x)  exp 


(- 
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4 I vST 


x -/5Sq  eX?  C-  ^ 

Taking  logarithms  on  Equation  (3-14)  yields 


(3-22) 


ln{  — exp  (tt^—  ) } 
mp  ^ v 2mp ' 


1 + u 
1 - u 


£L_ 

mp 

mp 


- 

mp 


(3-23) 


Cancelling  u in  Equations  (3-22)  and  (3-23)  gives 

1 


PVx)  - Wi  - exp  {-  In2 
m 2 q x r 2q 


mp  exp  ( - 


2mp' 


(3-24) 

It  can  be  seen  that  mp  exp  (-  q/2mp)  is  the  teometrical  mean 
of  x in  the  above  equation.  To  treat  the  distribution  of 
transmitted  irradiances  passing  through  a turbulent  medium 
when  the  nature  of  photons  is  also  considered,  Equations 
(2-60)  and  (3-24)  should  be  used  in  the  Mandel  formula . 

From  Equation  (3-24),  when  m is  very  large,  the  mean  and 
standard  deviation,  respectively,  are  given  by 

E(x)  = mp  (3-25) 

and 


"2 

°x  = nip { exp  (g-)  - 1}  (3-26) 

When  m is  very  large, 

ax  = 2mpq  = /E(x)q  (3-27) 

When  p is  very  small,  the  signal- to-noise  ratio,  S/N,  is 
given  by 
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S 

N 


(3-28) 


When  p is  small,  and  m is  extremely  large,  Equation  (3-24) 


Thus,  the  four  distributions,  Bernoulli's,  the  normal, 
Poissonian  and  the  lognormal  are  related.  The  lognormal 
distribution  has  the  advantage  of  Bernoulli's  and  Poissonian 
distributions  and  that  of  the  normal  distribution,  i.e.,  the 


In  the  previous  chapter  (Chapter  II) , the  signal-to- 
noise  is  independent  of  the  source  irradiance  in  a random 
medium.  In  this  chapter,  the  nature  of  photons  has  been 
taken  into  consideration  in  a homogeneous  medium,  and  the 
result  derived  is  that  the  signal- to-noise  ratio  increases 
as  the  source  irradiance  increases. 


becomes 


P (x) 

mv  ' 


(3-29) 


skewedness  and  the  convenience  in  calculation,  respectively. 


CHAPTER  IV 

THE  LOGNORMAL  EMISSION  PROFILE 


Gaussian  Doppler  Profile 


140-145 


According  to  the  principles  of  special  relativity,  when 
a light  source  is  approaching  an  observer,  the  expression 
for  the  Doppler  shift  of  frequency  is  given  by 


where  v is  the  observed  frequency  of  the  radiation  from  a 
light  source  which  approaches  the  observer  with  a velocity 
along  the  direction  of  the  sight,  v , v is  the  observed 
frequency  of  the  radiation  of  the  light  source  when  it  is  at 
rest,  and  c is  the  velocity  of  light  for  a medium  with  refrac- 
tive index  of  unity. 

When  v <<  c,  Equation  (4-1)  becomes 


v 

v 


o 


(4-1) 


V 

V 


(1  + — ) 
c ' 


V 

X 


(4-2) 


o 


or 


X 

X 


v 


(1  - 


c 


-) 


(4-3) 


o 


where  X and  Aq  are  the  observed  wavelengths  of  the  radiating 
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source  in  motion  and  at  rest,  respectively.  If  the  light 
source  is  receding  from  the  observer,  one  merely  changes  the 
sign  of  the  velocity,  v , in  Equation  (4-1),  (4-2),  or  (4-3). 
The  expression  for  the  Doppler  profile  can  be  obtained  from 
Equations  (4-2)  and  (4-3) , and  the  Maxwellian  distribution 
of  the  kinetic  velocities  along  the  sight,  P(v  ),  given  by 


p(vx>  = (set)'2  exP  <-  2FT>  <4-4> 

where  m is  the  mass  of  the  emitting  particle,  k is  the 
Boltzmann  constant,  and  T is  the  absolute  temperature  of 
the  medium. Therefore,  the  expression  for  the  Doppler 
profile,  P(v)  or  P(A),  on  the  scale  of  frequency  or  wavelength, 
respectively,  is  given  by 


mv 


dv 


P(v)  = P(vx)|^| 


( 


2 ! 2 
me  N % r me 

-j)  exp  { - 


2TTkTv 


o 


2kT 


(v 


,2, 

V } 


o 


(4-5) 


or 


P(A)  = P(v) 


dv 


x i 


x'  1 dA 


= ( 


2 2 
me  \ % r me  , , . ,2, 

7)  exp  { 7 (X  - A ) } 

2^kTX“  2kT  Z ° 


(4-6) 

Equation  (4-5)  or  (4-6)  shows  that  the  Doppler  emission 
profile  is  symmetric  about  the  mean.  However,  careful  ex- 
perimental studies  of  the  emission  profiles  showed  that  they 

1/q 

were  not  symmetrical.  Integrating  Equations  (4-5)  and 

(4-6)  for  v and  X from  -°°  to  °°,  respectively,  gives 

CO 

/.„P(v)  dv  = 


1 


(4-7) 
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and 

CO 

i.ja)  dx  = 1 (4-3) 

Equation  (4-7)  or  (4-8)  indicates  that  the  probability 
of  photons,  emitted  from  the  sources  with  the  thermal  motions, 
distributed  in  the  range  for  v or  A from  -°°  to  °°,  respective- 
ly, is  unity.  But  "negative  wave  numbers"  (frequencies  or 
wavelengths)  "have  no  physical  reality . "d44  if  the  integration 
is  done  in  the  interval  from  0 to  «>  for  v or  A , in  accordance 
with  the  physical  reality,  then  the  probability  is  less  than 
unity  which  is  against  the  physical  reality.  Hence,  the 
average  energy  of  emitted  photons,  E , calculated  from 
Equation  (4-5)  or  (4-6) , is  given  by 

oo 

Eave  = /0h  dv  > hvo  (4-9) 

or 


Eave  “ hor  P<X>  dX  * ” > O'lO) 

It  is  interesting  to  note  that  (i)  the  average  energy  of  the 
photons  emitted  from  sources  with  thermal  motions  is  greater 
than  the  average  energy  of  photons  emitted  from  sources  at 
rest,  hvQ  or  hc/AQ,  as  shown  in  Equation  (4-9)  or  (4-10) 
respectively,  and  (ii)  the  emission  profile  in  terms  of  the 
wavelength  as  shown  in  Equation  (4-10)  is  unreasonable. 

The  half-widths  of  the  Doppler  broadenings  derived  from 
Equations  (4-5)  and  (4-6)  are  given  byd4”* 


2v 


Av 


D 


o ,2kT  In  2s  % 


c 


m 


(4-11) 
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and 


AX 


D 


2Ao  /2kT  In  2v% 
c ' m ' 


(4-12) 


Lognormal  Hypothesis 

The  hypothesis  is  based  upon  the  following  reasons: 

(i)  In  spectroscopy,  both  frequency  and  wavelength  are 
positive  quantities.  Similarly,  variables  are  positive 

in  the  lognormal  distribution  as  shown  in  Equation  (2-35). 

(ii)  In  the  emission  profile,  the  width  of  a spectral  line 
is  proportional  to  the  mean  value  (frequency  or  wave- 
length) as  shown  in  Equation  (4-11)  or  (4-12).  Simi- 
larly, in  the  lognormal  distribution,  the  standard 
deviation  is  proportional  to  the  mean  as  shown  in 
Equation  (2-39).  In  the  normal  distribution,  the 

standard  deviation  and  the  mean  should  be  statistically 
1 46 

independent . 

(iii)  In  the  emission  profile,  the  expected  energy  of  a 
random  system  is  greater  than  the  expected  energy  of 
a homogeneous  system  as  shown  in  Equations  (4-9)  and 
(4-10).  Similarly,  in  the  lognormal  distribution,  the 
expected  value  of  a variable  in  a random  system  is 
greater  than  the  expected  value  of  a variable  in  a 
homogeneous  system,  as  shown  in  Equation  (2-38) . 

Hence,  it  would  be  reasonable  to  assume  that  the  Doppler 
profile  would  take  a lognormal  shape,  given  by 
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P(v)  = 


/2tt  3dv 


exp  (~ 


23 


D 


In2  — ) 
o 


(4-13) 


or 


P(A)  = 


-=T In2  f-) 

/2?  3dA  23^  o 


(4-14) 


where  3D  is  a dispersion  parameter  for  the  spectral  broadening 
due  to  the  thermal  randomness. 

The  probability  of  photons  emitted  from  sources  experienc- 
ing thermal  motions,  in  the  range  from  0 to  ®,  for  v or  A , 
is  given  by 


/ P(v)  dv  = 1 

o 


(4-15) 


or 


/ P (A)  dA  = 1 


(4-16) 


The  average  energy  of  photons  from  sources  experiencing  thermal 
motion,  Eave,  obtained  from  Equation  (4-13)  or  (4-14)  is 
given  by,  respectively, 


or 


00  3 

!ave  = /0h^p(v)  dv  = hvQ  exp  (^) 


co  g ^ 

Eave  = /0"TP(A)  dX  = T~  exp 


(4-17) 


o 


(4-18) 


It  is  seen  from  the  above  expressions,  that  the  average 
energy  of  photons  from  the  sources  experiencing  thermal 
motions  is  greater  than  that  of  photons  from  the  sources  at 
rest.  A discussion  will  be  given  to  understand  where  the 
excess  energy  comes  from  in  the  above  expressions  and  to 
evaluate  Bp. 
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We  can  consider  that  half  of  the  emitting  particles  are 
moving  toward  the  observer  and  the  other  half  of  the  emitting 
particles  are  receding  from  the  observer  with  the  root-mean- 
square  velocities  v^  and  -v  , respectively,  which  are  given 
by 


1-2  1 
2mvx  = IkT 


(4-19) 


The  average  energy  of  photons  given  off  by  the  emitting 
particles,  Eave>  from  Equation  (4-1),  is  given  by 


ave 


= hv 


hv  [ (- 
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1 + — i 
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-)*] 
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-2 

- hvo  d - 
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(4-20) 


The  above  equation  indicates  that  the  average  energy  of 

photons  from  the  sources  experiencing  thermal  motions,  E 

° ’ ave ’ 

is  greater  than  that  of  photons  from  the  sources  at  rest  hv 

’ o 

Substituting  Equation  (4-20)  into  Equation  (4-17)  or 
(4-18), 


( D*  ,, 

exp  (^~)  = (1 


-2 

v i 

_Xs  -% 
2> 


(4-21) 


and  taking  logarithms  on  both  sides  of  the  above  equation 
yields 


-2 

2 v 

% = ~ In  [1  - -f] 

cz 


(4-22) 


For  vx  <<  c ’ in  tke  analytical  laboratory,  then  the  above 


equation  becomes 
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X\  2 


- ^ 


(4-23) 

Substituting  Equation  (4-19)  into  the  above  equation  gives 


kT 


me 


(4-24) 


Substituting  Equation  (4-24)  into  Equations  (4-13)  and  (4-14) 
gives 


P(v)  = (— 

27TkTv 


me2 

'WT 


7)  2 exP  In''  ^— ] 


.2  v_ 
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(4-25) 
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and 


P(A)  ~ ( — mC  2)^  exP 


2ukTA' 


me  , 2 A , 

~2kT  in  A~J 
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(4-26) 


The  above  equations  show  that  the  Doppler  emission  profiles 
are  not  symmetric  and  the  integrations  for  v and  A from  0 to 
°°  are  unity.  It  has  been  shown  that  the  lognormal  distribution, 
Equation  (2-35),  is  approximated  by  the  normal  distribution, 

r\ 

Equation  (2-42),  for  small  values  of  3 . Similarly,  the  log- 
normal Doppler  emission  profile,  Equation  (4-25)  or  (4-26), 
can  be  approximated  by  the  normal  Doppler  emission  profile, 
Equation  (4-5)  or  (4-6),  for  small  values  of  kT/mc2.  The 
asymmetry  and  the  peakedness  of  the  lognormal  Doppler  emission 
profile,  Equation  (4-25)  or  (4-26),  can  be  easily  calculated 
from  Equations  (2-40)  and  (2-21),  respectively. 

So  far,  the  advantages  of  the  lognormal  expressions  over 
the  Gaussian  expressions  have  been  shown.  The  mathematical 
validity  of  the  lognormal  expression  will  be  given  in  the 
^^xt  chapter  to  show  that  the  Equations  (4-25)  and  (4-26)  are 
more  accurate  than  Equations  (4-5)  and  (4-6) . 


Lorentzian  Profile 


Another  major  cause  of  line  broadening  in  a gas  results 
from  collisions  of  radiating  particles  (atoms,  ions  or 
molecules)  with  each  other.  The  following  expression  describe 
the  frequency  distribution  of  collision-broadened  (Lorenthian) 
line 


P(v)  - (§£) 


(v  - Vo)2  + (|P2 


(4-27) 


where  vq  is  the  center  frequency,  and  Av  is  the  width  between 


the  half-power  points  of  the  curve 
The  integration 


147 


J P(v)  dv  < / P(v)  dv  = 1 


(4-28) 


indicates  that  the  probability  of  a photon  distributed  in 
the  range,  v from  0 to  is  less  than  unity.  The  average 
energy  emitted  from  a Lorentzian  line,  given  by 

oo 

Eave  = /0hv  P<v)  dv  = oo  (4-29) 

is  not  well  defined,  because  the  integration  is  divergent . 
Beside  the  conflicts  to  the  physical  reality,  the  Lorentzian 
line  cannot  explain  the  shift  and  asymmetry  of  the  line  pro- 
file. It  can  only  describe  that  the  line  profile  is  more 
peaked  than  the  Gaussian  profile. 
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Lognormal  Shape  of  Doppler-Collisional  Broadening 


Assuming  that  the  Doppler 
broadening  are  two  independent 
Doppler-collisional  broadening, 

P(v)  = — - — exp  ( — Ky 

/2¥3v  2gz 


broadening  and  the  collisional 
processes,  the  profile  for  the 
P (v) , is  given  by 

ln  — ) (4-30) 

os 


where 


= e 


D 


+ 


(4-31) 


and  Bq  is  the  dispersion  parameter  of  the  spectral  profile 
due  to  the  collisional  process,  and  vQg  is  the  geometrical 
mean  of  the  frequencies. 

Assuming  the  energy  of  the  system  of  interest  remains 
unchanged  due  to  the  collisional  process,  the  average  energy, 
emitted  from  a Doppler-collisional  profile,  is  given  by 


E 


ave 


/ hv  P(v)  dv 
0 


- hvos  exp  + 1 g2} 

Equations  (4-17)  and  (4-32)  give 


(4-32) 


vos  vo  exP  (4-33) 

Similarly,  the  Doppler-collisional  profile  in  terms  of 
wavelength  is  given  by 


P(A)  - — - — exp  ( — kj-  In 

/ O n ~\  r\ 


2 A 


/2tt  3A 


23' 


os 


(4-34) 


os  the  geometrical  mean  of  the  wavelengths, 


where  A 


given  by 


Aos  ~ Ao  exP  3C)  (4-35) 

Because  the  reliable  interpretations  of  the  shifts  in 
the  spectral  profiles  were  mostly  done  with  the  absorption 
method,  hence  the  shifts  of  the  spectral  profiles  will  be 
discussed  in  Chapter  VI. 


CHAPTER  V 

THE  DOPPLER  EFFECT  AND  THE  DISTRIBUTION 
OF  VELOCITIES 

In  this  chapter,  an  approximate  expression  for  the 
Doppler  effect  is  discussed.  The  lognormal  distribution  of 
emission  profile  is  also  derived  from  the  approximate  expression 
of  the  Doppler  effect  with  mathematical  validity.  The  dis- 
tributions of  thermal  velocities  are  also  easily  derived  from 
the  lognormal  distribution  of  emission  profile  without  the 
concept  of  phase  space  being  involved. 


Approximate  Expression  for  Doppler  Effect 


It  is  assumed  that  the  approximate  expression  for  the 
Doppler  effect  is  given  by 


In  ±- 
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(5-1) 


or 
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exp  (^) 


(5-2) 


When  vx  « c,  expanding  the  above  equation  yields 
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(5-3) 
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The  polynomial  expansion  of  Equation  (4-1) , derived 
from  the  principles  of  relativity,  when  v <<  c,  is 

X 

.2 

(5-4) 


v 


V 1 v 

v 1 + c + J 

o c 


It  is  seen  from  Equations  (5-3)  and  (5-4)  that  the  Equation 
(5-1)  is  more  accurate  than  the  Equation  (4-2) 


y_ 

v 


= 1 + 


v 


Equation  (5-1)  and  the  Maxwell  distribution  of  one-dimensional 
velocitied  of  Equation  (4-4)  give 


P(v)  = 


1 / 1 

exP  ( 7 

/2tFBv  23 


In2  — ) 
v J 
o 


(5-5) 


where  3 = (kT/mc2)^. 

The  above  equation  is  identical  with  Equation  (4-25) , 
which  is  derived  from  the  lognormal  assumption. 

It  is  concluded  that  the  lognormal  expression  for  the 
Doppler  profile,  Equation  (4-25),  is  more  valid  than  Gaussian 
expression,  Equation  (4-13),  because  Equation  (5-1)  suffers 
less  approximation  than  Equation  (4-2). 


The  Distribution  of  Velocities 


In  the  studies  of  plasmas,  the  flux  of  particles 
reaching  the  electrodes  has  been  represented  by  three  different 
expressions  for  the  same  thermal  velocity  in  treating  the  same 
problem. ^51  Considerations  should  be  based  upon  the  dis- 


62 


tributions  of  velocities.  The  mathematical  derivations  for 
the  distributions  of  velocities  are  quite  complicated.  How- 
ever they  can  be  easily  derived  from  the  lognormal  approach. 

The  one-dimensional  distribution  of  velocities,  P(v  ), 
can  be  easily  derived,  from  the  lognormal  Doppler  profile,  in 
Equation  (4-25),  and  the  Doppler  effect,  in  Equation  (5-1), 
and  is  given  by 

P<V  - pcv)  ig^-i 

x 

2 

z m N%  , mvxN 

- hiFh  exP  (-ZkT5  (5-6) 

The  two-dimensional  distribution  of  velocities,  which  is 
defined  as  two  one-dimensional  distributions  of  velocities 
in  x and  y directions,  sometimes  called  the  circular  normal 

of  Rayleigh  distribution,  can  be  easily  derived  as 

„ 2 
vo  v? 

P(v2)  = ~ t exP  ( 7)  (5-7) 

8 8 

where  8 - (kT/mc  ) 2,  = (vx  4-  Vy) 2 and  v is  the  velocity 

of  the  particle  along  the  y-axis . 13 ’ 152 

The  three-dimensional  distribution  of  velocities  in 
x,  y and  z directions,  sometimes  called  the  sperical  normal 
or  Maxwell  distribution,  can  be  easily  derived  as 

p(v  ) = £L  (_|)  exp  ( — |)  (5_8) 

/i  8Z  282 

where  8 . = (kT/mc  )2,  v-,  = (v2  + v2  + v2)^  and  v is  the 

^ y £ z 

velocity  of  the  particle  along  the  z-axis . 32 > 152 


CHAPTER  VI 

LOGNORMAL  ABSORPTION  PROFILES 


Introduction 


In  an  actual  absorption  experiment,  the  intensity  of 
light  transmitted  through  a fixed  layer  of  material  is  re- 
corded while  the  frequency  of  the  incident  light  is  varied. 
From  the  assumption  that  the  decrease  of  the  light  intensity, 
I,  on  passage  through  material  of  thickness,  dx,  is  propor- 
tional to  I dx,  it  is  inferred  that  the  variation  of  I with 
depth  x follows  the  law 

Vx>  = v°>  exP  (_kvx)  (6-1) 

where  Iv ( 0)  is  the  incident  intensity  at  frequency  v (Hz)  and 

kv  ^cm  "*")  is  caH-ed  absorption  coefficient  at  frequency  v.1^7 

In  order  to  evaluate  quantitatively,  it  is  advisable 

to  take  a look  at  the  process  of  absorption.1^7  Consider  a 

parallel  beam  of  light  of  frequency  between  v and  v + dv  and 

intensity  1^  travelling  in  the  positive  x-direction  through 

a layer  of  atoms  bounded  by  the  planes  x and  x + dx . For 

simplicity,  let  the  refractive  index  of  this  medium  be  unity, 

and  the  velocity  of  light  be  c.  The  phase  front  will  travel 

through  the  slab  of  thickness  dx  in  the  time  dt  = dx/c. 

Suppose  there  are  N,  atoms/cm3  in  level  1,  of  which  dN,  are 

Iv 
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capable  of  absorbing  the  incident  radiation  and  dN2v  excited 
atoms  capable  of  emitting  the  same  frequency  radiation.  Neg- 
lecting the  effect  of  spontaneous  re-emission  in  view  of  the 
fact  that  it  occurs  equally  in  all  directions,  the  decrease 
in  energy  of  the  incident  beam  is  given  by 


-d(Ivdv)  = hv(B12dNlv 


B21dN2v>1 


dx 

VC 


(6-2) 


where  B-^  is  the  coefficient  of  stimulated  absorption  for  the 
transition  from  level  1 to  level  2 and  B^  is  the  coefficient 
of  stimulated  emission  for  the  transition  from  level  2 to 
level  1. 

From  the  above  equation,  it  follows  that 


1_ 

I 

v 


dl 

v 

dx 


dv 


hv 

c 


(B12dNlv 


B21dN2v> 


(6-3) 


Recognizing  that  the  left-hand  member  is  k^dv  as  defined  by 
Equation  (6-1) , the  above  equation  becomes 

kvdv  = ir  <B12dNlv  - B21dN2v>  <6-4> 


Assuming  that  coefficients  B-^  and  B^  are  independent 
of  frequency,  then  integrating  over  the  entire  absorption 
line  yields  approximately 


/ k dv 
J o v 


/ k dv 

J „ v 


(b12Ni 


b21n2) 


(6-5) 


where  vQ  is  the  frequency  at  the  center  of  the  line,  assuming 

that  the  atoms  are  normally  distributed  about  v . 

o 


When  the 
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stimulated  emission  is  ignored,  then  the  above  equation 
becomes 


hkvdv  ■ “E2  B12N1  (6-6) 

A better  integration  of  Equation  (6-4) , assuming  that 
the  stimulated  coefficients  are  independent  of  frequency  and 
that  the  stimulated  emission  is  negligible,  is  given  byd^d 


I 


o v 


— dv  = - B,  0N, 

c 12  1 


For  the  atomic  line,  Equation  (6-6) 

154 


of  Equation  (6-7) . 
tion  (6-7)  is  favorable 


However  for 
155 


(6-7) 

is  a good  approximation 
the  molecular  band,  Equa- 


Lognormal  Absorption  Profile 


Assuming  that  the  distribution  of  photons  follows  the 
lognormal  model,  as  shown  by  Equation  (4-30),  we  have 

dNlv  = N1  p(v)dv  (6-8) 

Assuming  that  the  stimulated  emission  is  negligible,  and 
substituting  the  above  equation  into  Equation  (6-4)  yields 


ir  Bi2Nip(')) 


Using  the  relation 

.2 


B 


12 


ire 

m hv 
e 


(6-9) 


where  e and  mg  are  the  charge  and  the  mass  of  the  electron, 
respectively,  and  f is  the  oscillator  strength  between  level 
1 and  level  2,  then  Equation  (6-9)  becomes 
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"e  N,  f — exp  (~\  In2  Z— ) 

26z  v 


v in  c 1 /t% — _ 

e /2tt3v 


(6-11) 


os 


Integrating  the  above  equation,  assuming  f is  constant, 
because  oscillator  strength  is  dimensionless,  yields 

2 


/ k dv  = N,  f 

J o v m c 1 

e 


(6-12) 


Defining  vmax  as  the  frequency  where  is  maximum,  as 


dk 

v 

dv 


= 0 


(6-13) 


From  Equation  (4-33) , the  above  equation  gives 


= v exp  (-3; 


3 r2\ 
2 


(6-14) 


'max  vo  ' PD 

The  above  equation  indicates  that  both  the  Doppler  and 
collisional  effects  will  cause  a red-shift  of  the  absorption 
line.  Ignoring  the  slit  function  of  a grating  monochromator 
when  the  slit  is  very  narrow,  then 


VA 


= k 


ne  N f ^ \ 


m c 
e 


) exp  ( — — In2  y — ) 
/2?3c  23  os 


Integrating  the  above  equation  yields 

2 


/ k dA 
0 A 


m c' 
e 


N!  f Xo  exP  <7  6D  + 


(6-15) 


(6-16) 


The  above  equation  indicates  that  the  absorption  integral 
(or  integrated  absorption  coefficient)  increases  with  in- 
creasing 3q  and  3^ • Defining  ^max  as  the  wavelength  where 
k^  is  maximum  as 
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0 


(6-17) 


then 


A 


max 


(6-18) 


The  above  equation  indicates  that  both  the  Doppler  and 
collisional  effects  will  cause  a red-shift  of  an  absorption 
line.  Because  Equation  (6-4)  is  valid  for  molecular  bands, 
assuming  3^  increases  with  increasing  temperature,  one  will 
expect  the  experimental  result  showing  a red-shift  of  molec- 
ular absorption  bands  in  the  gaseous  phase,  when  temperature 

1 C4 

becomes  higher. 

Both  the  temperature  and  the  collision  will  affect  the 
emission  and  the  absorption  profiles,  as  shown  in  Equations 
(6-11)  and  (6-15).  Hence,  for  example,  in  the  atomic  fluo- 
rescence spectroscopy  with  a line  source,  the  fluorescence 
signal  or  the  quantum  yield  of  an  excitation  process  will  be 
different  if  the  temperature  is  changed,  or  if  the  collisional 
effect  is  changed,  even  if  there  is  not  collisional  excitation 
or  de-excitation. 


It  will  be  shown  here  that  the  lognormal  absorption 
profiles  can  be  approximated  by  the  Lorentzian  absorption 
profiles . 

The  Lorentzian  absorption  profiles,  with  frequency  and 
wavelength  scales,  respectively,  are  given  by 


Lognormal  and  Lorentzian  Absorption  Profiles 
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k = k 
v o 


Av 

2 7 

4(v  - v ) + Av 


(6-19) 


and 


k,  = k 
A o 


AA 

4(A  - *"7  + A A 2 


(6-20) 


where  kQ  (cm  ) is  the  value  of  kv  (cm-1)  or  k^  (cm-1)  at 
v = vq  or  A = Aq , respectively,  and  Av  or  AA  is  the  half- 
width of  the  Lorentzian  absorption  profile,  Equation  (6-19) 
or  (6-20),  respectively.157 

Ignoring  the  small  effect  of  spectral  shift  in  Equation 
(6-11)  gives 

2 


kv  - <5^>  N1  f <-^-> 

e 


/21t3v  exp  In2  — ) 

23Z  Vo 


When  v is  in  the  vicinity  of  v , the  above  equation  becomes 


ku  = Ni f <— — ) 


/2tt8v  1 + — y In2  — 
° o o 2 v 

23  o 


= N1  f ^-4 ) 

e /2tt  3v 


26  vo 

Because  — kQ  at  v = v^,  the  above  eauation  can  be  expressed 


as 


k = k 
v o 


2 2 
83  vz 
o 


4 (v  - vQ)2  + 882v2 


Let  Av  - 2/28vo,  then  the  above  equation,  which  is  derived 
from  the  lognormal  absorption  profile,  Equation  (6-11),  be- 


69 


comes  the  Lorentzian  absorption  profile,  Equation  (6-19) 
when  v is  in  the  vicinity  of  v . When  v becomes  far  from  v 

U Q > 

Equations  (6-11)  and  (6-19)  approach  zero.  Hence,  the  log- 
normal absorption  profile,  Equation  (6-11),  is  approximated 
by  the  Lotentzian  profile,  Equation  (6-19),  when  written  in 
terms  of  frequency. 

Similarly,  when  written  in  terms  of  wavelength,  the 
lognormal  absorption  profile,  Equation  (6-15),  can  be  approx- 
imated by  the  Lorentzian  absorption  profile,  Equation  (6-20). 

The  Lorentzian  absorption  profiles  written  in  terms  of 
frequency  and  wavelength  units,  Equations  (6-19)  and  (6-20), 
respectively,  are  symmetric  about  their  maxima.  It  can  be 
seen  that  if  an  absorption  profile  is  symmetric  on  the 
frequency  scale,  then  the  absorption  profile  reolotted  versus 
wavelength  cannot  be  symmetric,  because  of  the  reciprocal 
relationship  between  frequency  and  wavelength.  Therefore, 
the  Lorentzian  absorption  profiles  with  frequency  and  wave- 
length scales,  Equations  (6-19)  and  (6-20),  cannot  both  be 
mathematically  valid. 


CHAPTER  VII 

ERRORS  IN  THE  REACTIONS  OF  FIRST-ORDER 


Error  of  the  Average  Life-Time 

Suppose  we  have  a reaction  of  first-order  given  by 

A * x + " (7-1) 

then,  the  reaction  rate  can  be  expressed  as 


d[A] 

dt 


k [A] 


(7-2) 


where  [A]  is  the  concentration  of  the  species  A (cm-3,  M 
or  torr)  at  time  t (s) , and  k is  the  first-order  rate  constant 
of  the  reaction  (s-'*')  . 

Integration  of  the  above  equation  gives 


[A]  = [A] ^ exp  (-kt) 


(7-3) 


where  [A] ^ is  the  initial  concentration  of  species  A. 


The  average  life-time,  t,  and  the  standard  deviation  of 


the  life-time  of  species  A,  a , can  be  expressed  as 


T = 


uu 

ht(nr-> 

dt 


158,159 


(7-4) 


and 


la 

= {-° 


x)2(^I)  dt 


1^ 
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/.  <^>  ^ 


(7-5) 
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Using  the  formula 

oo 

.n  _ n! 


/ tn  exp  (-st)  dt 


(7-6) 


and  substituting  Equation  (7-6)  into  Equations  (7-4)  and 
(7-5),  yields159’160 


T = k 


(7-7) 


1 

k 


T K (7-8) 

From  the  above  equations,  the  relative  errors  of  the  life- 
time in  a first-order  reaction  is  given  by 

~ - 100%  (7-9) 

The  above  equation  is  used  to  evaluate  the  natural 
broadening  of  a spectral  line  where  the  life-time  of  an  ex- 
cited atom  is  equal  to  the  standard  deviation  of  the  life- 
time.15^ Because  no  reliable  measurements  of  life-time  were 
reported  with  relative  errors  of  100%,  it  is  possible  that 
the  above  equation  is  only  valid  for  a system  consisting  of 
one  particle. 

For  a system  consisting  of  m particles,  the  relative 
error  of  the  life-time  may  be  approximately  given  by 


x,m  _ 100  , 


% 


(7-10) 


m vdn 

where  x is  the  average  life-time  of  n particles  and  o is 
111  x ,m 

the  standard  deviation  of  life-time  in  a system  of  n particles 
For  example,  in  order  to  have  a relative  error  less  than  1%, 
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ignoring  other  instrumental  errors,  the  system  of  interest 
should  consist  of  tenthousand  particles  as  a minimum. 

Time  Dependent  Rate  Constant 


A study  in  a real  example  will  be  considered,  as  in  the 
thermal  decomposition  of  azoisoprane  at  270°C,  given  by161 
(CH3 ) 2^H=NCH (CH3 ) 2 -*■  N2  + ^6^14  (7-11) 

Assuming  that  the  reaction  is  first-order,  the  initial 
pressure  of  azoisoprane  is  (torr) , and  at  time  t,  its 
pressure  becomes  pA  (torr),  with  a decrease  of  x (torr),  then161 
PA  = Pi  - x (7-12) 

From  the  rate  equation 


dp, 

dt 


-kp 


A 


(7-13) 


then 


k = r In 


Pi  - x 


(7-14) 


In  the  above  equation,  the  partial  pressure  of  azoiso- 
prane is  not  measurable,  only  the  total  pressure,  p,  is 
measurable,  which  is  the  summation  of  the  three  partial 
pressures  given  by 

P = Pa  + Pn2  + pc6h14  = Pi  + x (7-15) 

Cancelling  x in  Equations  (7-14)  and  (7-15)  gives 


k - \ In  Pl 

t 2p.  - p 


(7-16) 
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The  experimental  results  for  the  reaction  in  Equation  (7-11) 
are  duplicated  in  Table  7-1.  ^ It  is  interesting  to  see 
that  the  calculated  rate  constant  decreases  with  increasing 
time.  In  order  to  explain  this,  a new  model  is  proposed 

this  model,  it  is  assumed  that,  due  to  some  unknown 
parameters,  e.g.,  the  microscopic  inhomogeneity  of  the  tem- 
perature, the  distribution  of  the  rate  constant,  P(k),  follows 
the  normal  law,  given  by 


P(k)  = — — exp  [ — ^ (k  - k )2] 

/Z?B  2?  ° 


(7-17) 


where  0 is  the  dispersion  parameter  (s  2)  of  the  rate  constant 
due  to  the  microscopic  inhomogeneity  of  the  rate  constant 
and  kQ  is  the  theoretical  rate  constant  or  the  macroscopic 
average  of  the  rate  constants. 

If  the  rate  constant  is  microscopically  homogeneous , the 
theoretical  total  pressure,  p at  time,  t,  is  given  by 


k - i In  Pi 
0 t 2pi  - p 


o 


(7-18) 


From  Equations  (7-16)  and  (7-17),  the  probability  func- 
tion  of  the  total  pressure  at  time  t,  P(p) , is  given  by 


p(p>  = paoi^l 


1 

/2rr0t(2pi  - p) 


exp 


r 1 i 2 

l Tj—rj  In 

2 ZtZ 


(7-19) 

Therefore,  the  macroscopic  (measured)  total  pressure  at  time 
t,  p,  is  given  by 
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Table  7-1 

The  Thermal  Decomposition  of  Azoisopropane 

t (s) 

P (torr)  k x 103  (s-1)161 

0 

35.15  

180 

46.30  2.12 

360 

53.90  2.11 

540 

58.85  2.07 

720 

62.20  2.03 

1020 

65.55  1.96 

average  2.02 
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— oo 


2 2 

= 2Pi  - (2Pi  - pQ)  exp  (^-|-) 


(7-20) 


The  calculated  rate  constant,  k,  is  then  given  by 


The  above  equation  shows  that  the  calculated  rate  constant 
decreases  with  increasing  time,  as  shown  in  Table  7-1.  A 
least-square  fitting  of  the  above  equation  based  on  the  data 
in  Table  7-1  is  shown  in  Figure  7-1.  The  time  dependent 
rate  constant  is  given  by 


Both  the  experimental  result  and  the  theory  confirm 
that  the  calculated  rate  constant  is  not  a constant,  but  a 
time  dependent  variable. 


k 


2.17  x 10 


1.99  x 10"7  t 


(7-22) 


Measured  Rate  Constant 
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Fig.  7-1  - Time  Dependence  of  the  Rate  Constant 


CHAPTER  VIII 

GENERAL  DISCUSSIONS  ON  THE 
LOGNORMAL  DISTRIBUTION 


Normal  (Gaussian)  Distribution 


Suppose  the  variable  z is  distributed  according  to  the 
normal  distribution  given  by 

i 2 

P(z)  = exp  (-|-)  , for  -oo  < z < oo  (8-1) 

/2tt 

"It  follows  that  the  normal  model  (distribution)  is  strictly 
applicable  only  if-  • -negative  measurements  are  possible. 
"But  many  experimental  variables  are  always  positive  (e.g., 
volume,  mass,  concentration,  kinetic  energy,  absolute  tem- 
perature, and  the  entropy  change  of  an  isolated  system).  In 
observations  of  variable  of  this  sort,  the  Gaussian  distribu- 
tion function  is  often  inappropriately  applied  - the  experi- 
mental distribution  is  usually  positively  skewed."16  The 
discussion  for  the  distribution  of  the  positive  variables  is 
given  in  the  following  section. 

1 7 

Two-Parameter  Lognormal  Distribution1 

The  two-parameter  lognormal  distribution,  sometimes 
simply  called  lognormal  distribution,  is  written  as 
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P(x)  = < 


— - — exp  ( — In  ^— ) , for  0 < x < 
/2?0x  2B2  xo 


(3-2) 

L 0 , for  x <.  0 
where  xq  > 0 and  xQ  is  the  geometrical  mean  of  the  variable 
x.5  The  above  equation  was  obtained  from  Equation  (8-1) 
through  the  transformation 


z = i In  *- 
6 xo 


(8-3) 


In  this  distribution,  one  of  the  most  important  properties  is 
that  the  standard  deviation  is  proportional  to  the  mean  of 
the  variable  x as  shown  in  Equation  (2-39) . It  is  interesting 
that  many  positive  physical  quantities  can  be  described  by 
this  property.  For  example,  the  error  in  the  weighing  of 
a sample  is  usually  assumed  to  be  proportional  to  the  weight 
in  the  analytical  laboratory.  In  other  words,  the  heavier  the 
sample,  the  greater  the  error.  In  the  quantitative  analysis 
laboratory  of  the  Department  of  Chemistry,  the  errors  of 
students'  reports  are  assumed  to  be  proportional  to  the  true 
values  of  their  unknown  samples.  Therefore,  the  grading  of 
students'  reports  of  the  unknown  samples  is  based  upon  their 
relative  errors.  In  analytical  chemistry,  the  precision  is 
usually  expressed  as  the  ratio  of  the  standard  deviation  to 
the  mean,  i.e.,  the  proportionality  constant  between  the 
standard  deviation  and  the  mean,  e.g.,  the  relative  error,  or, 
as  the  ratio  of  the  mean  to  the  standard  deviation,  e.g.,  the 
signal- to-no is e ratio  or  the  spectral  resolving  power. 

Another  property  in  this  distribution  is  that  the  logarithm 
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of  a positive  variable  follows  the  normal  distribution.  Many 
positive  physical  quantities  are  expressed  in  the  logarithmic 
forms,  e.g.,  volume,  pressure,  concentration,  activity  coef- 
ficient in  thermodynamics,  wavelength,  frequency  and  distance 
between  electron  and  the  nucleus.162"166  The  distribution 
of  particle  sizes  is  one  of  the  best  applications  of  the 
lognormal  distribution.  The  distribution  of  particle  sizes 
is  usually  given  by 


exp  ( — Kj  In2  ~)  , for  0 < r < 


P (r)  = ^ /^3rr  2Br  "o 

ta0,  for  r < 0 (8-4) 

where  Br  is  a dispersion  parameter  of  the  particle  radius,  r 
is  the  radius  of  the  particle  (m) , and  r^  is  the  geometrical 


mean  of  the  radii  (m) , or  the  expected  radius  where  Br  = 0 . 
The  above  equation  was  derived  from  Equation  (8-1)  through 
the  transformation6 


33 


z = 


h ln  - 

6r 


(8-5) 


The  lognormal  distributions  for  the  cross  sections,  Q 

2 

(m  ),  and  the  weights,  W (kg),  of  particles  can  be  derived 
from  Equation  (8-4)  and  written  as 


P(Q)  = 


/2Tp(2Br)Q 


exp  [- 


1 2 Q , 

7 ln  0~] 

2(2B  ) 'O 


and 


(8-6) 


P(W)  = 


/2tF  (3  Br)W 


exp  [ In2  H_] 

2(3B^)  Wo 


(8-7) 
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through  the  transformations 
Q = Trr2 


Q = irr 
o o 


TT  4 3 

W = j -rrr  p 

TT  4 3 

Wo  " J 1Irop 

where  p is  the  density  of  the  particle  (kg/nr) 
Define 
3q  = 26r 


(8-8) 

(8-9) 

(8-10) 

(8-11) 


= 3 3 

W r 


(8-12) 

(8-13) 


then 


P(Q)  = 


P(W)  = 


/2¥3qQ 


/2?3WW 


exp  ( — ln2  o~) 
2e2  Qo 


exp  (—^r  In2  ~) 


23 


(8-14) 


(8-15) 


'W 


From  Equations  (8-4),  (8-14)  and  (8-15),  it  is  seen  that  if 

the  particle  radii  follow  the  lognormal  distribution,  then 

the  cross  sections  and  the  weights  should  also  follow  the 

lognormal  distribution.  However,  it  was  shown  that  if  the 

particle  radii  followed  the  normal  distribution,  then  the 

cross  sections  or  the  weights  would  not  follow  the  normal 
5 21 

distribution.  ’ From  Equations  (2-40)  and  (2-41),  it  is 
known  that  if  the  dispersion  parameter  is  small,  the  normal 
distribution  is  an  approximation  of  the  lognormal  distribu- 


tion . 


18 


Because  3r  is  smaller  than  3y,  and  if  the  normal 
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approximation  is  to  be  used,  then  the  normal  distribution 
gives  a better  approximation  for  the  particle  sizes  than  for 
particle  weights.  Yuan  used  the  normal  distribution  for 
students'  heights  and  the  lognormal  distribution  for  their 
weights  in  his  theses.^  This  can  be  easily  justified  if  it 
is  assumed  that  the  weight  is  roughly  proportional  to  the 
third  power  of  the  height.  From  Equations  (8-5),  (8-6)  and 
(8-7)  one  can  find  the  following  expressions  given  by 


Q ^ irr 
ave  ave 


and 


W ^ Trr^  p 

ave  3 ave 

where  ? ave , Qave  and  W&ve  are  the  average  radius,  cross 
section  and  weight  of  the  lognormally  distributed  particles. 

It  is  interesting  to  know  that  the  above  expressions  do  hold 
in  the  manipulation  of  data,  as  shown  in  the  following  example 
Suppose  there  are  two  spherical  particles  with  radii  1 and  3 
(cm).  Therefore,  their  cross  sections  are  tt  and  qir  (cm2), 
respectively.  The  average  radius  and  cross  section  are, 
respectively,  given  by 


rave  2 


Qave  = 5 17 


(8-16) 


O 

It  is  seen  from  the  above  equations  that  Q it r 

n xave  ave 

The  concentrations  of  many  trace  elements  in  geological 
samples  were  found  to  follow  the  lognormal  distribution  given 
by 
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p(c>  = — = — exP  (— K ln2  £“)  (8-17) 

/2?3C  2BZ  Co 

where  C is  the  concentration  (M) , and  6 is  the  dispersion 

7 o 77 

parameter  of  the  concentration. 

If  the  concentration  of  excited  atoms  follows  the  log- 
normal distribution,  then  the  emission  intensity  in  atomic 
emission  spectroscopy,  should  follow  the  lognormal  distribu- 
tion. Equation  (8-17)  can  be  understood  if  it  is  assumed  that 
the  distribution  of  free  energies  of  the  samples  follows  the 
normal  distribution.  If  it  is  assumed  that  the  concentration 
of  a material  of  interest  (atom,  ion  or  molecule)  in  an 
optical  mdeium  follows  the  lognormal  distribution  given  by 
Equation  (8-17),  the  distribution  of  the  transmitted  irradiances , 
P(I),  through  this  medium,  can  be  expressed  as 


P(I) 


1 

/??BI  In  (I00/I) 


exp 


T 

~oo 
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(8-18) 

through  the  relations 


ln 

(Ioo^  2.303e  1 C 

(8-19) 

ln 

(Ioo/V  2 . 303e  1 Co 

(8-20) 

where  e is  the  molar  absorptivity  (M_1cm_1)  and  1 is  the 
length  of  the  optical  medium  (cm) , Iqo  is  the  incident 
irradiance  and  I is  the  transmitted  irradiance  when  the 
medium  is  homogeneous . 

If  3 in  Equation  (8-17)  is  very  small,  then  Equation 


(8-17)  becomes 
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P(G)  = — — 1 — exp  [ i-j(C  - C )2]  (8-21) 
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and  the  distribution  of  the  transmitted  irradiances  follows 
the  two-parameter  lognormal  distribution  given  by 


P(I)  = 
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The  above  equation  is  identical  with  Equation  (2-35)  except 
for  the  notations  for  the  dispersion  parameter. 

The  lognormal  distribution  can  also  be  applied  to  the 
study  of  the  propagation  of  indeterminate  error.  In  the 
case 


R = AB/C  (8-23) 

one  would  inquire  to  "how  scatter  in  measurements  of  quantities 
A,  B,  C etc. , is  translated  into  random  variation  in  the 
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final  result  R."  In  the  normal  distribution,  it  is  assumed 

that  the  errors  are  quadraticallv  additive  for  independent 
processes.  But  in  this  case,  the  error  of  R is  expressed  as 
a. 


<T>2  - + 4b2  + <!r>2  + 


(8-24) 


where  a denotes  the  standard  deviation,  assuming  that  the 
standard  deviations  are  approximately  equal  to  the  square 
roots  of  the  variances . 

It  is  interesting  to  note  that  in  Equation  (8-24) , the 
assumption  is  that  the  relative  errors,  instead  of  the  absolute 
errors,  are  quadratitally  additive.  If  it  is  assumed  that 
the  lognormal  distribution  is  applicable  to  Equation  (8-23) , 
i.e.,  the  logarithms  of  the  quantities  follow  the  normal 


distribution,  then  we  have 


In  R — In  A + In  B + In  C + 


(8-25) 


and 
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2 , 2 
In  A + aln  B 


(8-26) 


assuming  that  the  errors  of  the  logarithms  of  the  quantities 
are  quadratically  additive.  If  the  errors  are  small,  Equation 
(8-24)  is  found  from  Equation  (8-26). 

To  evaluate  the  precision  of  an  instrument,  it  is  con- 
venient to  express  the  noise  of  a variable  in  terms  of  its 
logarithmic  form.  For  example,  let  us  evaluate  the  noise 

caused  by  a recorder.  The  signal,  H,  is  assumed  to  be  given 
by 


where  H is  the  final  output  signal  (cm) , R is  a conversion 
factor  of  the  recorder  (cm/V) , and  V is  the  output  voltage 
of  a detector  (V) . 

The  above  equation  gives 


H 


R V 


(8-27) 


aH  = <Vctr)2  + (R°v)2  % 


(8-28) 


If  the  noise  from  the  recorder,  N^,  is  expressed  as 
NR  * V°R 


(8-29) 


and  the  overall  noise  from  the  detector,  Ny,  as 


(8-30) 


then  the  final  output  noise,  , is  given  by 


a 


H 


K + N2)% 


(8-31) 


The  above  equation  is  used  in  the  analytical  laboratory, 
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assuming  that  NR  and  are  independent . However,  it  is  seen 
from  Equations  (8-29)  and  (8-30)  that  NR  and  Ny  are  not  inde- 
pendent at  all.  But,  if  we  assume  that  the  logarithms  of 
the  quantities  in  Equation  (8-27)  follow  the  normal  distribution, 
so  from  Equation  (8-27) 

In  H = In  R + In  V 

then 

0 In  H = aln  R + aln  V (8-32) 

then  a R and  a y are  the  independent  noises.  Hence, 
o in  h s^ould  be  used  to  characterize  the  precision  of  the 
recorder  function. 

Three-Parameter  Lognormal  Distribution 


Recently,  a three-parameter  lognormal  distribution, 

P(I'),  for  the  irradiance  of  an  optical  wave  passing  through 

the  turbulent  medium  was  proposed  and  experimentally  veri- 
123 

fied.  The  equation  is  given  by 
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for  I < I- 


(8-33) 

where  I'  is  the  measured  irradiances , It  and  I are  both 

1 o 

positive  constants,  Iq  is  the  geometrical  mean  of  (I  - 1-^)  , 
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and  1^  is  called  the  lower  limit  of  the  transmitted  irradiances 
This  distribution  is  also  called  the  lognormal  distribution 
with  lower  limit.  It  was  derived  empirically  from  Equation 
(8-1)  through  the  transformation5 

I'  - I-, 
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In  (- 


0 


o 


(8-34) 


Equation  (8-33)  can  also  be  derived  from  Equation  (2-35)  by 
the  relation 

I'  = I + Ix 

where  I'  is  the  measured  irradiance,  I is  the  transmitted 
irradiance  of  the  light  beam  from  the  source,  and  1-^  is  the 
background  irradiance  which  is  a constant. 

In  a gas  chromatograph,  the  retention  time  of  a sample, 
t,  should  have  a lower  limit,  t which  is  the  retention 

a. 

time  of  the  air  peak.  Thus,  the  shape  of  a sample  peak  in 
the  gas  chromatograph  might  be  reasonably  expressed  as 


P(t) 


— 1 exp  [-1  In2  (— -■  -ta)l  (8-35) 
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where  tQ  is  the  geometrical  mean  of  (t  - t ) and  g is  a dis- 
persion parameter  of  the  sample  for  the  broadening  of  a peak. 
The  above  equation  is  plotted  in  Figure  8-1.  The  curves  in 
Figure  8-1,  generated  from  a simple  empirical  equation,  are 
similar  to  those  shown  in  Figure  8-2  proposed  by  Li.168  A 
theoretical  curve  based  on  the  binomial  (Bernoulli's)  dis- 
tribution is  duplicated  in  Figure  8-3,  for  comparison,  showing 

that  the  shape  of  a sample  peak  is  skewed  and  has  a lower 
169 

limit.  It  is  seen  that  the  three-parameter  model  is  good 


and  convenient. 


. 15 
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Fig.  8-1  - Shape  of  the  Gas  Chromatographic  Peaks  (t  = 100;  t = 126) 
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u=  3 cm 
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Fig.  8-2  - Shape  of  the  Gas  Chromatographic  Peaks 


Fraction  of  Solute 
Leaving  Column 


Number  of  Plate  Gas  Volumes 
Leaving  Column 


Fig.  8-3  - Elution  Curve  for  One  Unit  of  Solute 
From  Five-Plate  Column 


The  time,  t 

’ max’ 

(8-35),  is  given  by 


where  the  neak  is  maximum  in  Equation 
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t = t + t exp  (-B; 

max  a o K v / 


(8-36) 


and  the  maximum  peak  height,  h,  is  therefore 

, „2 


h = 


/2tt6  t 


exp  (£-) 


(8-37) 


The  area  under  the  curve,  A,  is  unity.  Using  the  expression 
for  the  number  of  theoretical  plates,  N,  given  by 
ht 


N = 2ir  (-■  ™ax)2 


(8-38) 


and  from  Equations  (8-36),  (8-37)  and  (8-38),  the  number  of 
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theoretical  plates  is  given  by 


N = exp  (82)  t + t exo  (-32)  (8-39) 

3ztz  a ° 

o 

The  above  equation  is  plotted  in  Figure  8-4.  It  is  shown 
in  Figure  8-4  that  the  number  of  theoretical  plates  decreases 
then  increases  with  increasing  3. 


Four-Parameter  Lognormal  Distribution 


The  four-parameter  lognormal  distribution  is  an  extension 
of  the  normal  distribution  to  allow  for  both  a lower  and  upper 
limit  to  the  possible  values  of  the  variable,  given  by 
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Fig.  8-4  - Number  of  Theoretical  Plates. 
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for  x,  < x < x 
1 - - u 

P(x)  = 0,  for  x ^ x or  x < x 

1 u 

(8-40) 

where  x^  and  x-^  is  the  upper  and  the  lower  limit  of  the 
variable  x,  and  (xQ  - xq)/(xu  - x ) is  the  geometrical  mean 
of  (x  - x1)/ (xu  - x) . 

The  above  equation  was  derived  from  Equation  (8-1) 
through  the  transformation 

-i  x - x,  x - x, 

z * (ln  s— ' ln  sjpr-jb  (8‘41) 

Equation  (8-40)  was  used  to  study  the  distribution  of 
the  percentages  of  the  lawns  covered  by  clover,  by  letting 
x denote  the  percentage  of  coverage,  xQ  = 0 , and  xu  = 100. 

It  was  used  to  study  the  distribution  of  the  droplet  sizes 
with  an  upper  limit  of  droplet  size,  by  letting  x denote  the 
diameter  of  the  particle  and  x-^  = 0 . The  above  empirical 
approach  was  experimentally  shown  to  be  better  than  the  em- 
pirical Rosin-Rawmler  equation,  and  the  semi-empirical  Nuki- 
yama-Tanaswa  equation . ^ > 171 

Folded  Lognormal  Distribution 

This  distribution  was  derived  from  the  model  of  random 
walk  with  reflecting  barrier  for  the  distribution  of  the  drop- 
let size  with  an  upper  limit,  given  by 
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(8-42) 


where  x is  the  droplet  diameter  (m)  and  x is  the  maximum 
diameter  (m).2^>132 

If  the  droplet  size  decreases  due  to  the  evaporation, 
the  droplet  diameter  at  any  time,  xfc , is  given  by 


x2  = x2  - kt  (8-43) 

where  k (m2  s ^)  is  an  evaporation  constant.^-22 

The  distribution  of  the  droplet  sizes  at  any  time  will 
be  given  by 
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(8-44) 

where  the  maximum  diameter  x^  ^ and  the  geometrical  mean  of 


diameters,  x 


at  any  time  are  given  by,  respectively 
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Equation  (8-44)  can  be  used  for  the  time  dependent  distribution 
of  droplet  sizes  in  the  atomic  flame  spectroscopy. 

Relationships  among  the  Lognormal  Distributions 

A three-parameter  lognormal  distribution,  as  in  Equations 
(8-33)  or  (8-35),  can  be  reduced  to  a two-parameter  log- 
normal distribution,  as  in  Equation  (8-4),  if  the  lower  limit 
is  zero,  i.e.,  1^  = 0 or  t = 0 . A four-parameter  lognormal 
distribution,  such  as  in  Equation  (8-40),  can  be  reduced  to 
a three-parameter  lognormal  distribution,  as  in  Equation 
(8-33)  or  (8-34),  if  the  upper  limit  is  infinite  and  it  can 
further  be  reduced  to  a two-parameter  distribution  if  the 
upper  limit  is  infinite  and  the  lower  limit  is  zero.  Similar- 
ly, a folded  lognormal  distribution,  such  as  in  Equation 
(8-42) , can  be  reduced  to  a two-parameter  lognormal  distribu- 
tion if  the  upper  limit  is  infinite. 


CHAPTER  IX 

INTERFERENCE  ON  MAGNESIUM  BY  TRACE  CONCOMITANTS 
IN  FLAME  ATOMIC  ABSORPTION  SPECTROMETRY 

Concomitant  interferences  remain  one  of  the  major  sources 
of  inaccuracy  in  flame  atomic  absorption  spectrometry.  Most 
previous  efforts  emphasize  interferences  due  to  large  excesses 
of  concomitants  in  the  sample  matrices.  In  this  chapter,  the 
interferences  by  the  presence  of  trace  concomitants  (on  mag- 
nesium) , which  have  been  generally  ignored  in  flame  atomic 
absorption  spectrometry  are  reported. 

Experimental 

Observations  were  made  with  Perkin-Elmer  Model  303 
atomic  absorption  spectrometer  (Perkin-Elmer,  Norwalk,  Conn.) 
with  a 10  cm  burner  for  an  air/acetylene  flame.  A Ca-Mg-Al 
hollow  cathode  lamp  (Varian  Techtron,  Palo  Alto,  Calif,) 
was  used  in  all  studies;  the  image  of  the  light  source  was 
3 mm  in  diameter  on  the  focal  plane;  the  monochromator  was 

O 

adjusted  to  2852  A for  magnesium  with  1 mm  slit  width.  A 
magnesium  stock  solution  was  prepared  directly  by  dissolving 
reagent  grade  magnesium  sulfate;  an  aluminum  stock  solution 
was  prepared  from  reagent  grade  aluminum  chloride  and 
gravimetrically  determined.  Serial  dilutions  were  carried 
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with  deionized  water.  In  several  studies,  an  electric  field 
at  50  V cm  was  applied  across  the  sprayer  and  the  burner 
head  of  the  nebulizer-burner. 


Results  and  Discussion 


Analytical  Growth  Curves 

The  analytical  curves  for  magnesium  are  depicted  in 
Figure  9-1  for  several  flame  heights  and  for  acetylene  and 
air  flow  rates  of  4.3  and  28.7  1 min-1,  respectively.  The 
slopes  of  the  analytical  curves  at  1.5  mm,  2 . 5 mm  and  5 mm 
are  less  than  unity.  In  the  low  region  (1.5  - 5.0  mm)  of 
the  flame,  magnesium  does  not  have  enough  time  travelling 
upward  in  the  flame  to  receive  sufficient  thermal-collisional 
energy  to  atomize.  The  slow  rate  of  volatization  is  most 
likely  responsible  for  the  departure  of  the  slopes  from  unity. 
According  to  the  Nukiyama-Tanasawa  expression,  the  concentra- 
tion of  MgSO^  in  solution  does  not  affect  the  size  of  sprayed 
droplets  in  the  nebulizing  chamber.111  The  lowering  of  the 
vapor  pressure  of  the  droplet  due  to  the  presence  of  salt  in 
the  droplet  can  be  qualitatively  understood  by 


r r .2cm  N ■,  r , 
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(9-1) 


where  Pr  is  the  equilibrium  vapor  pressure  of  water  over  a 

_ o 

droplet  of  radius  r (cm)  and  density  p (g  cm  ) containing 
m (g)  of  the  dissolved  MgSO^,  P^  is  the  equilibrium  vapor 
pressure  over  a flat  surface  of  water  at  the  same  temperature 
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Concentration  of  MgSO^  (lO-1^  mole/1) 


Fig.  9-1  - Analytical  Curves  of  Magnesium 
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T (°K),  (dyne  cm  ) is  the  surface  tension  of  water  at 

temperature  T,  I is  van't  Hoff's  factor,  W and  M are  the 

molecular  weights  of  MgSO^  and  water,  respectively,  and  R is 

173 

the  gas  constant.  From  the  above  thermodynamic  equation, 

it  can  be  seen  that  the  vapor  pressure  of  water  over  the 
droplet  can  be  reduced  significantly  by  the  presence  of 
dissolved  salt,  if  the  size  of  the  droplet  is  extremely  small. 
Therefore,  the  rate  of  desolvation  would  be  expected  to  be 
significantly  reduced  by  the  presence  of  salt  in  the  droplet, 
when  the  droplet  does  not  have  enough  time  to  evaporate  the 
solvent  in  the  low,  cool  region  of  the  flame.  After  the 
desolvation  is  completed,  the  size  of  dried  magnesium  par- 
ticulate will  depend  upon  the  original  solution  concentration. 
After  desolvation,  the  resulting  particulate  melts.  Large 
particulates  should  result  in  a lesser  lowering  of  the  melting 
point  than  small  particulates  according  to  the  Thomson  equa- 
tion 


2aV 

Lr 


(9-2) 


where  Tm  is  the  lowering  of  melting  point  of  a spherical 

solid  particulate  relative  to  the  macroscopic  melting  point 

T , a (dyne  cm  ) and  V (cm  mole  ) are  the  surface  tension 

and  molar  volume  of  the  dried  particulate  at  the  melting 

point,  respectively,  and  L (erg  mole  is  the  heat  of  melt- 
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mg. 


The  decrease  in  vapor  pressure  of  a large  molten  droplet 


is  expressed  by  the  Kelvin  equation 
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where  the  physical  parameters  are  those  of  the  molten  drop- 
let. Also,  a large  molten  droplet  results  in  slow  vapor- 

ization because  of  the  small  surface/mass  ratio.  All  of  the 
successive  factors,  low  vapor  pressure  of  solution  (water) 
droplets,  high  melting  point  of  solid  particulates,  low 
vapor  pressure  and  low  surface/mass  ratio  of  molten  drop- 
lets for  a droplet  formed  initially  from  a solution  of  high 
concentration,  are  present  and  cause  the  reduction  in  the 
slope  of  the  analytical  curve  of  magnesium  in  the  low,  cool 
region  of  the  flame.  The  slope  of  the  analytical  curve  at 
a 10  mm  measurement  height  (above  burner  top)  is  about  unity. 
The  slope  of  the  analytical  curve  at  a 20  mm  height  is  less 
than  unity  but  larger  than  those  at  1.5  mm  and  2 . 5 mm  measure- 
ment heights. 

At  a measurement  height  of  20  mm,  magnesium  salt  particles 
have  more  time  travelling  upward  in  the  flame  to  allow  for 
the  processes  of  solvent  evaporation  and  solute  vaporization. 
The  rate  of  volatization  of  magnesium  salts  would  appear  to 
be  relatively  unimportant.  The  negative  deviation  of  the 
analytical  curve  from  a slope  of  one  is  probably  a result  of 
the  lognormal  distribution  of  the  transmitted  intensities 
due  to  the  turbulence  of  the  flame  at  20  mm.  If  it  is  assumed 
that  the  transmitted  intensities  of  a beam  of  radiation  can 
be  described  by  Equation  (8-22) , then  the  average  transmitted 
intensity  measured,  I , is  given  by 

cL  V 0 
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where  I is  the  transmitted  intensity  when  there  is  no  tur- 
bulence of  the  concentration  of  magnesium  atoms  in  the  optical 
region  of  the  flame,  e is  the  molar  absorptivity  (M  ^ cm  ^) 
of  magnesium  atoms,  1 (cm)  is  the  horizontal  length  of  the 
flame  at  20  mm,  3 is  a dimensionless  parameter  for  the  dis- 
persion of  transmitted  intensities  and  Cq  (M)  is  the  concen- 
tration of  magnesium  atoms  when  the  flame  is  homogeneous. 

From  the  above  equation,  the  average  (measured)  absor- 
bance of  the  inhomogeneous  region  of  the  flame  is  given  by 
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where  I is  the  incident  intensity, 
oo  J 

It  is  seen  in  the  above  equation  that  the  measured 
absorbance  will  show  a negative  deviation,  due  to  the  tur- 
bulence of  the  flame,  as  C increases. 

’ o 


Interferences 

Interferences  of  magnesium  in  solutions  of  10  ^ M MgSO^ 
and  various  concentrations  of  AlCl^  are  shown  in  Figure  9-2. 
Flow  rates  are  the  same  as  those  used  for  the  data  in  Figure 
9-1.  The  interference  curve  of  Mg  by  AlCl^  is  generally  flat 
at  a 20  mm  height  above  the  burner  top.  The  matrix  effect 
is  visible  at  high  concentrations  of  A1C1~.  It  is  reasonable 
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Fig.  9-2  - Interference  Curves  of  Magnesium  by  Aluminum 
Chloride 


to  conclude  from  the  curve  at  20  mm  that  the  presence  of 
AlCl^  for  the  range  of  Mg  concentrations  studied  does  not 
affect  the  size  of  droplets  sprayed  and  the  nebulization 
efficienty  of  the  solutions . ^ The  interference  curves  at 
1.5  mm  and  2.5  mm  look  peculiar  but  are  in  good  agreement 
with  that  obtained  by  Menzies.^^  However,  the  present 
curves  are  more  detailed  in  the  low  concentration  range  of 
AlCl^.  It  seems  obvious  that  the  peculiar  shapes  of  inter- 
ference curves  in  the  lower  region  of  the  flame  is  due  to 
the  presence  of  some  unknown  eutectic  mixtures  of  magnesium- 
aluminum  compounds  and  the  mass-controlled  processes  of 
vaporization.  Measurements  at  5 mm  and  10  mm  heights  result 
in  higher  signals  than  at  20  mm  but  suffer  more  interference 
Practically,  it  would  be  advisable  to  gain  accuracy  at  the 
sacrifice  of  signal  intensity  by  looking  at  higher  parts  of 
the  flame. 

The  data  in  Figure  9-3  were  obtained  by  changing  the 
flow  rate  of  acetylene  and  by  applying  the  electric  field 
across  the  sprayer  and  burner  head.  The  flow  rate  of  air 
was  the  same  as  mentioned  above.  Higher  flow  rates  of  acety 
lene  enhanced  the  signal  at  a 20  mm  height  but  reduced  the 
signal  at  a 2.5  mm  height  when  no  electric  field  was  applied 
Higher  flow  rates  of  acetylene  caused  the  transit  time  of 
magnesium  to  be  shorter.  Therefore,  signals  are  lower  and 
interferences  are  greater  at  lower  heights  with  higher  fuel 
flow  rates.  Also,  higher  flow  rates  of  acetylene  would  re- 
sult in  a larger,  more  reducing  flame  and  so  the  absorption 
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L 

L 

H 


Fig.  9-3  - Effects  of  Acetylene  Flow  Rate,  Flame  Height 

and  Electric  Field  on  the  Interference  Curves 
H:  Acetylene,  51/min.;  L:  Acetylene,  4.3  1/min. ; 
-o- : 20  mm  above  Burner  Head;  -x- : 2.5  mm  above 
Burner  Head;  N:  Normal  Operation;  F:  With  Applied 
Electric  Field. 
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signals  are  larger  and  interferences  are  less  in  higher 
parts  of  the  flame.  To  reduce  the  interference,  the  slightly 
fuel-rich  flame  is  favorable  with  measurements  in  the  higher 
parts  of  the  flame.  The  parallelism  between  the  interference 
curves  with  and  without  electric  field  indicates  that  the 
decreased  signals  when  the  electric  field  was  applied  might 
be  due  to  the  increased  rates  of  coalescence  and/or  deposition 
of  the  sprayed  droplets  on  the  walls  and  spoiler  of  the 
nebulizing  chamber  by  means  of  electrostatic  ef f ects . 

The  interference  curves  of  magnesium  by  AlCl^  in  the 
presence  of  0.1  M HC1  were  generally  flat  and  slightly  bent 
down  at  high  concentration  of  AlCl^ . To  reduce  chemical 
interferences  due  to  trace  concentration  of  AlCl^,  HC1  should 
be  added  to  the  measured  solutions . 


DETERMINATE  ERRORS  IN  UV-VISIBLE 
MOLECULAR  ABSORPTION  SPECTROSCOPY 


CHAPTER  X 

DETERMINATE  ERRORS  IN  ULTRAVIOLET  AND  VISIBLE 
SPECTROPHOTOMETRIC  MEASUREMENTS  CAUSED  BY  MULTIPLE 
REFLECTIONS  WITHIN  THE  CELL 

Absorption  spectrophotometry  has  been  a widely  used 

analytical  tool.  Departures  from  Beer's  law  due  to  multiple 

reflections  within  the  cell  have  been  Quantitatively  approxi- 

1 78-1 8 2 

mated  in  the  past  by  other  workers.  The  departures 

were  expressed  in  terms  of  absorbance  errors  and  correction 

1 7 8-1 82 

factors  for  the  measured  fraction  of  light  absorbed. 

Positive  deviation  of  absorbance  should  result  in  a correction 

factor  of  less  than  unity  for  the  measured  fraction  of  liyht 

absorbed.  However,  contradictions  arise  when  attempting  to 

correlate  the  positive  absorbance  errors  and  correction  fac- 

178-1 82 

tors  of  greater  than  unity.  To  remove  the  confusion 

and  to  allow  a better  understanding  of  the  factors  leading 
to  greater  accuracy  in  absorption  spectrophotometry,  a more 
rigorous  approach  to  the  effects  caused  by  multiple  reflections 
within  the  cell  is  required. 

Theoretical  Considerations 

When  radiation  passes  through  an  interface  between  two 
transparent  semi-infinite  media  of  different  refractive  indi- 


106 


107 


ces,  a fraction  of  the  radiation  is  reflected  at  the  inter- 
face. If  the  incident  angle  is  zero,  as  is  normally  the  case 
in  absorption  spectrophotometry,  the  fraction  of  light  re- 
flected is 


f 


nl  " n2 
nl  + n2 


(10-1) 


where  n^  and  are  the  refractive  indices  of  the  two  media 
comprising  the  interface.  When  a parallel  beam  of  monochro- 
matic radiation  with  zero  incident  angle  passes  from  one 
transparent  semi-infinite  medium  to  another  transparent  semi- 
infinite medium  through  a window  of  internal  transmittance, 
t , the  actual  fractions  fo  the  light  transmitted,  reflected 
and  absorbed  due  to  multiple  reflections  of  the  radiation 
within  the  parallel  surfaces  of  the  window  are  t-^.  r^  and 
a^2>  respectively.  In  Figure  10-1  the  fractions  transmitted, 
reflected  and  absorbed  are  shown.  By  considering  the  processes 
of  transmission,  reflection  and  absorption  at  the  two  inter- 
faces, the  fractions  transmitted,  t^>  reflected,  r-^2'  an<^ 
absorbed,  a-^.  are  given  by 
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12 
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12 


= - fi>a  - V'1  + hh<  + •••) 

(1  - gm  + f2tw)(i  - tw) 

1 - hVw 


(10-4a) 

(10-4b) 


Similarly,  when  the  radiation  passes  in  the  opposite 
direction  from  medium  2 to  medium  1 through  the  same  window 
with  zero  incident  angle,  the  fractions  of  radiation  trans- 
mitted, t2l’  reflected,  r2p>  an<^  absorbed,  a2p>  are  also  given 
by 

(1  - f1)(l  - f2)t; 


'21 


12 


w 


1 - f,f„tA 

1 2 w 

f,f2  - 2fnf0t2  + f0 
1 2 1 2 w 2 

1 - f -|  f q t 2 
1 2 w 


a 


12 


a + W*1  - f2)(1  - t«> 
1 - flVw 


(10-2b) 


(10-5) 


(10-6) 


The  above  equations  show  that  only  transmittance  is  symmetric 
with  respect  to  the  direction  of  light. 

In  order  to  develop  realistic  expressions  for  absorption 
spectrophotometry  consider  an  absorption  cell  consisting  of 
two  identical  parallel  windows  separated  by  a parallel  column 
of  absorbing  medium  (the  sample  which  is  assumed  to  be  homo- 
geneous) of  transmittance,  t . When  a parallel  beam  of  mono- 
chromatic radiation  passes  through  the  cell  with  zero  incident 
angles,  the  fractions  of  light  transmitted,  T , reflected,  Rt , 
and  absorbed,  at , due  to  multiple  reflections  within  the  cell 
(see  Figure  10-2)  are  given  by 
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1 - r21ts 

Rt 

= r12  + r21t12t21ts(^1  + r21ts  + ■'')  (10-8a) 

Rc 
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21  s 

To  simplify  the  derivations,  the  following  assumptions 
are  made: 


(i) 

183 

the  radiation  is  monochromatic; 

(ii) 

the  cell  windows  are  flat,  have  parallel  sides,  and 
have  no  imperfections; 

(iii) 

the  optical  alignment  allows  only  zero  incident 
angles  at  the  interfaces  of  interest; 

(iv) 

the  temperature,  pressure  and  reflactive  index  are 
the  same  in  the  light  path  and  in  the  surrounding 
of  the  path;  thermo-optic  effects  are  absent; 

(v) 

The  absorber  is  isotropic,  non-f luorescent , non- 
black-body  radiative  and  homogeneous;  the  absorber 
results  in  no  Bragg  diffraction,  no  dichroism  and 

112 

no  Mie  scattering;  i.e,,  there  are  no  solid  or 


liquid  crystals  formed,  and  all  particles  are  much 

smaller  than  the  wavelengths  of  interest  in  the 

, , 184 

absorber ; 

(Vi) 

the  Fabry-Perot  interferometric  line  spacings  are 
too  small  to  be  resolved  by  the  detector ; 

(vii) 

Beer's  law  is  valid  at  all  concentrations;  there 
are  no  association,  dissociation,  and  cross-inter- 
action processes  which  cause  deviations;  the  ab- 
sorptivity coefficient  is  also  constant; 

(viii) 

when  two  cells  are  in  use,  they  are  identical; 

(ix) 

the  radiation  is  low  enough  to  guarantee  no  satura- 
tion of  excited  levels  of  the  absorbing  solute  and 

185 

there  are  no  other  systematic  errors; 

(x) 

finally,  window  absorption  of  commercial  fused 
silica,  solvent  absorption,  band  overlapping  of 
various  components  at  the  wavelength  of  interest, 
and  the  variation  of  refractive  index  of  the  solu- 
tion at  low  concentrations  are  all  assumed 

, . 186,187 

absent. 

Therefore,  Equations  (10-1)  through  (10-6)  can  be  reduced  to 


t = 

(1  - fT)(l  - f2) 

1 - f1f2  ’ (t  fc12  t21)  (1°  10) 

r 

fl  “ 2f lf 2 + f 2 

i - fp2  ; <r  = ri2  = r2i>  <10-n> 

where  t is  the  fraction  transmitted  for  one  of  the  cell  windows 

(t  = 1) , r is  the  fraction  reflected  for  one  of  the  cell 
w 
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windows  (t  = 1)  , and  a = 0 since  t =1  according:  to  the 
w w 

above  assumptions , and  therefore  Equations  (10-7)  through 
(10-9)  for  the  fractions  of  radiation  transmitted,  T , re- 
flected, R , and  absorbed,  a , when  radiation  passes  through 
the  cell  can  be  reduced  to 


R. 


A. 


i 2_2 

1 - r t 

s 


r + 


2 2 
rt  t 

s 

: T7L 

1 - r fcs 


(10-12) 


(10-13) 


a. 


(1  - r)(l  - tg) 
1 - rt 


(10-14) 


Absorbance  Error  and  Correction  Factors 


The  relative  absorbance  error  is  defined  as 


6A  = 


meas 


sol 


sol 


where  A , is  the  true  solution  absorbance  and  A 

sol  meas 


(10-15) 
is  the 


apparent  (measured)  absorbance.  Because  t = 1 for  the 
reference  cell  filled  with  pure  solvent,  the  apparent  (mea- 
sured) transmittance  of  the  solution  is  given  by 

<TAol 


'meas 


'meas 


(Tf) 


t' solv 


(1  - r )t; 

: 

1 - r Cs 


(10-16) 


where  (t,)  is  the  fraction  of  radiation  transmitted  through 
N t/sol 
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the  sample  cell  containing  sample  solution  and  (Tt)soqv  is 
the  fraction  of  radiation  transmitted  through  the  same  cell 
containing  solvent.  Hence,  the  apparent  (measured)  absorbance 
is  given  by 


A = A + log  ( 

meas  sol  & v 


i 2.2 

1 - r t 

s 


) 


(10-17) 


Equation  (10-15)  can  be  rewritten  as 

1 ( 1 ' r2  > 

l°g(- ^ 


r t 


6A  = 


log  (tg) 


(10-18) 


The  calculations  of  the  relative  absorbance  error  at 
25°C  were  carried  out  and  tabulated  in  Table  10-1  for  aqueous 

and  air  solutions  assuming  Corning  7940  fused  silica 

. , 188,189 

windows . 


Correction  Factor  for  Actual  Fraction 
of  Radiation  Absorbed 


An  expression  which  gives  correction  factors  for  the 
actual  fraction  of  light  absorbed  by  the  sample  is 


a. 


a. 


act 


a. 


meas 


1 - t 


(10-19) 


'meas 


i.e.  , 


act 


(1  - r) (1  + rt  ) 

TT7T— 


(10-20) 


is  de- 


where  the  measured  fraction  of  light  absorbed,  ameas > 
fined  as 

= 1 - t.___  (10-21) 


a 


meas 


meas 
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Table  10-1 

Relative  Absorbance  Error  (Parts  per  Thousand) 
Due  to  Multiple  Reflections 


Wavelength  (nm) 

275  468  656 


Absorbance  Aqueous  Air  Aqueous  Air  Aqueous  Air 
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Rewriting  Equation  (10-16) , we  have 


1 + /[(l  - r2)2  + 4r2t2  ] 

v ' meas 

o 2+. 

2r  t 


(10-22) 


meas 

The  calculations  of  the  correction  factor  for  actual  fraction 
of  light  absorbed  at  25°C  were  carried  out  and  tabulated  in 

Table  10-2  for  aqueous  and  air  solutions  for  fused  silica 

. , 188-190 

windows . 


Correction  Factor  for  Theoretical  Fraction 
~ of  Radiation  Absorbed 


Correction  factors  for  theoretical  fraction  of  light 
absorbed  by  the  sample  without  multiple  reflections  are  defined 
as 


1 - t 


theo 


a. 


meas 


1 - t 


meas 


and  by  use  of  the  previous  expressions  is  found  to 


theo 


i 2^2 

1 - r cs 

1 + A2 


(10-23) 


be 

(10-24) 


s 

The  calculations  of  correction  factors  for  the  theoretical 
fraction  of  light  absorbed  at  25°C  were  carried  out  and  tabu- 
lated in  Table  10-3  for  aqueous  and  air  solutions  for  fused 

...  . , 188,189 

silica  windows. 


Relative  Increase  of  Sensitivity  of  Absorbance 
Due  to  Multiple  Reflections 


The  relative  increase  of  sensitivity  of  absorbance  due  to 


Correction  Factors  for  Actual  Fraction  of  Light  Absorbed 
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multiple  reflections  within  the  cell  is  defined  as 


6 S 


dA 


"me  as 
dC 


dA 


sol 


dC 


dA 


sol 


dC 


(10-25) 


where  C is  the  concentration  of  the  absorbing  species  (the 
analyte)  in  solution  and  the  differentials  are  the  slopes  of 
the  analytical  curve  of  A or  A , versus  C,  respectively, 
and  A and  A were  previously  defined.  An  expression 

IT103.S  S O JL 

for  the  relative  increase  of  sensitivity  of  absorbance  is 
derived  as 


6 S = 


2r2t2 

s 

1 2 2 
1 - r t 

s 


(10-26) 


The  calculation  of  the  relative-  increase  of  sensitivity  of 
absorbance  at  25° C were  carried  out  and  are  tabulated  in  Table 

10-4  for  aqueous  and  air  solutions  for  fused  silica 

. , 188,189 

windows . 


Minimum  Concentration  Error 


The  concentration  error  caused  by  the  transmittance  error 
in  the  instrumental  system  due  to  multiple  reflections  is 
defined  as 


r - 1 rdc^ 

E " cTE 


(10-27) 


meas 

The  concentration  error  caused  by  the  transmittance  error  is 


derived  as 


Relative  Increase  of  Sensitivity  of  Absorbance  (Parts  per  Thousand) 
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E = 


(1  - r2t  )2 
s' 


TI 


(1  - r ) (1  + r t ) t In  (t  ) 

o S 


(10-28) 


The  minimum  values  of  concentration  error  (absorber  at 
25 °C)  are  tabulated  in  Table  10-5  for  aqueous  and  air  solutions 
for  fused  silica  windows  . ^8  ’ ^9 


Relative  Increase  of  Concentration  Error 
at  the  Same  Concentration 


The  relative  increase  of  concentration  error  at  the  same 
concentration  of  absorbing  species  caused  by  the  transmittance 
error  in  the  instrumental  system  due  to  multiple  reflections 
is  defined  here  as 


dt 


SE, 


meas 


(dC)  _ 1 

''C  ; dt 


(dC) 

'■C  ; 


dt 


(dC) 


(10-29) 


s 

An  expression  for  the  relative  increase  of  concentration  error 
at  the  same  concentration  is  then  derived  as 

(1  - r2t2) 2 

6Er  = « S ■»  o-  - 1 (10-30) 

L (1  - r)Z(l  + 4Ztp 

The  calculations  of  the  relative  increase  of  concentra- 
tion error  at  the  same  concentration  caused  by  transmittance 
error  in  the  instrumental  system  due  to  multiple  reflections 
(absorber  at  25 °C)  were  carried  out  and  are  tabulated  in  Table 

10-6  for  aqueous  and  air  solutions  for  fused  silica 
188,189 


windows . 


Minimum  Value  of  Concentration  Error 


121 


Relative  Increase  of  Concentration  Error  (Parts  per  Thousand) 

at  the  Same  Concentration 
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Relative  Increase  of  Concentration  Error 
at  the  Same  Transmittance 


The  relative  increase  of  concentration  error  at  the  same 
transmittance  caused  by  the  transmittance  error  in  the  instru- 
mental system  due  to  multiple  reflections  is  defined  as 


1 

HE  t 


6 E . 


meas 


Tn  (t  y 

meas  meas 


(10-31) 


t In  (t  I 

meas  meas 


An  expression  for  the  relative  increase  of  concentration  error 
at  the  same  transmittance  is  then  derived  as 


6E, 


i 2.2 

1 - r s 

1 + r2t2 
s 


tl  + A] 


(10-32) 


The  calculations  of  the  relative  increase  of  concentration 
error  at  the  same  transmittance  caused  by  the  transmittance 
error  in  the  instrumental  system  due  to  multiple  reflections 
(absorber  at  25° C)  were  carried  out  and  are  tabulated  in  Table 

10-7  for  aqueous  and  air  solutions  for  fused  silica 

. , 188,189 

mdows . ’ 


w 


Results  and  Discussion 


The  results  of  Table  10-1  are  in  good  agreement  with  those 
of  Burnett. Our  approach  is  easier  to  program  and  has  fewer 
simplifying  approximations.  The  results  in  Table  10-2  and 
Table  10-3  are  in  direct  conflict  with  the  results  obtained 
by  Hund  and  Hill,  Calvert  and  Pitts,  and  Jones  and  Salis- 
bury . Their  results,  which  we  feel  are  in  error, 


Relative  Increase  of  Concentration  Error  (Parts  per  Thousand) 

at  the  Same  Transmittance 
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have  resulted  from  their  rather  gross  approximations.  When 
windows  of  higher  indices  of  refraction  were  in  use,  lower 
limits  of  detection  (see  Table  10-1) , higher  sensitivities 
of  absorbance  (see  Table  10-4) , and  less  concentration  errors 
(see  Table  10-6)  could  be  achieved  due  to  the  multiple  reflec- 
tions in  absorption  spectroohotometry  for  trace  analysis. 
However,  at  higher  concentrations  of  absorbind  species  in 
solution,  concentration  errors  are  increased  (see  Table  10-6) , 
and  the  minimum  values  of  concentration  error,  when  multiple 
reflections  within  the  cell  are  present , are  higher  than  the 
conventional  minimum  value  of  concentration  error,  2.718, 
when  multiple  reflections  are  absent  (see  Table  10-5) . We 
feel  the  approach  we  used  can  be  anplied  to  the  study  of  the 

distortion  of  a molecular  absorption  band  due  to  the  frequency- 

184 

dependent  indices  of  water  and  of  the  cell  windows . 
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